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INTRODUCTION 


There  is  at  present  a  rapidly  expanding  body  of  empirical  data  and 
theoretical  Investigations  to  substantiate  the  concept  of  of  the  inter¬ 
action  between  the  geomagnetic  field  (GMF)  and  the  solar  corpuscular  radi¬ 
ation  (SCR)  and  its  concomitant  effects.  The  evidence  for  some  sort  of 
interaction  between  the  geomagnetic  field  and  the  tenuous,  highly  ionized 
plasma  issuing  from  the  sion  is  indisputable,  but  the  precise  nature  of 
this  interaction  is  still  an  open  question.  This  is  due  partly  to  the 
inherent  complexity  of  the  theoretical  problem,  and  pau-tly  to  the  lack  of 
sufficient  experimental  data.  However,  it  ceui  be  said  that  the  approxi¬ 
mate  theory  which  does  exist  agrees  generally  with  the  direct  measurements 
which  have  been  made  by  rocket  probes. 

The  general  outline  of  this  theory  is  that  the  solsu*  corona  is  con¬ 
tinually  expanding  outward  through  the  solar  system  in  the  form  of  a  very 
low-density  fully-ionized  plasma  consisting  of  protons  and  free  electrons. 
This  plasma  is  variously  known  as  the  interplanetary  gas,  the  solar  wind, 
and  the  solar  corpuscular  radiation  (SCR);  its  density  near  the  earth  is 
estimated  to  be  around  10  prbtons/cm;^  ,  drift  velocity  300  kn^sec  , 
and  thermal  kinetic  energies  between  10^  °K  and  10^  °K  .  It  is  gen¬ 
erally  believed,  without  any  notable  dissent  as  of  this  writing,  that  the 
GMF  is  almost  ccmpletely  separated  from  the  SCR  by  eun  interface,  or  trans- 
sition  region;  theoreticeLL  estimates  yield  a  minimum  thickness  of  the  order 
of  a  proton  Larmor  radius  in  the  ca4F  --  about  one  kilometer.  In  fact. 

Grad,  (1961),  has  shown  this  condition  to  obtain  generally  for  transition 
regions  between  magnetic  fields  and  pletsmas,  where  frozen- flow  donlnates 
over  the  diffusion  of  plasma  across  field  lines.  This  concept  of  the  In¬ 
terface  will  be  qualified  somewhat  as  we  proceed,  but  remains  essentially 
valid. 

The  Inqpact  of  the  SCR  congresses  the  GMF  on  the  day  side,  and  elongates 
It  on  the  night  side  of  the  earth.  Since  the  thermal  pressure,  nkT,  of 
the  SCR  Is  less  than  one-tenth  of  Its  dynamic  pressure  i  nmv^  ,  the  shape 
of  the  aC-SCR  Interface  Is  largely  Independent  of  the  temperature  of  the 
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stream  on  the  day  side.  However,  as  one  moves  away  frcxn  the  sun  the  GMF 
decreases  and. the  thermal  pressure  in  the  SCR  should  eventuslly  close  the 
interface,  so  the  field  lines  do  not  extend  to  infinity. 

The  interface  may  thus  he  considered  as  a  mathematical  surface  in  a 
free-houndary  problem  separating  field  on  one  side  of  the  boundary  from 
plasma  on  the  other.  The  boundary  conditions  then  require  surface  cur¬ 
rents,  J,  in  the  interface  which  must  exactly  cancel  the  (distorted)  GMF 
outside  of  the  surface  and  aijgment  it  inside  the  surface,  in  the  region 
known  as  the  magnetosphere.  If  there  is  to  be  no  field  outside  the  mag¬ 
netosphere  then  the  total  field  must  be  everywhere  tangential  to  the  sxu:- 
face. 

There  is  good  reason  to  believe  that  the  size  euid  the  shape  of  the 
interface  is  largely  determined  by  a  balance  of  forces  due  to  J  x  H 
force  on  the  surface  currents,  and  the  change  of  momentum  of  plasma  par¬ 
ticles  reflected  from  the  interface.  The  situation  is  depicted  schematically 
in  Fig.  1;  both  theory  and  measiirement  agree  on  r^  =  10  R  (R  =  earth 
radius ) 

In  addition  to  this  steady-state  effect  there  is  also  a  definite 
correlation  between  solar  activity,  geomagnetic  storms,  and  aviroral  activ¬ 
ity.  Particles  from  a  solar  flare  may  average  speeds  of  1000  km/sec  or 
more,  with  Increased  densities,  depending  upon  the  strength  of  the  flare 
(iGT,  1962a);  the  increased  flux  density  perturbs  and  ccmipresBes  the  GMF 
giving  rise  to  the  increeised  fluctuations  in  the  OtF  which  measiired 
at  the  earth's  surface.  Furthermore,  the  perturbation  of  field  lines  may 
cause  "dumping"  of  trapped  particles  from  the  radiation  belts  into  the 
auroral  zones,  which  gives  rise  to  increased  avuroral  displays  during  geo¬ 
magnetic  storms  (IGY,  196la). 

In  this  review  our  emphasis  will  be  primarily  on  the  theoretical 
aspects  of  the  problem,  especleOly  the  quantitative,  though  idealized, 
theory  which  has  been  worked  out  in  detail.  However,  there  is  also  a 
large  body  of  semi -quantitative  and  qualitative  theory  to  which  we  shall 
refer,  although  not  with  the  same  attention  to  mathematical  details.  Since 
it  is  impossible  to  review  all  the  materleCL  pertednlng  to  this  subject, 
our  aim  is  to  make  this  review  self-contained  by  presenting  enouc^  material 
on  each  aspect  to  provide  a  physlceG.  basis  and  theoretical  framework  for 
understanding.  , 
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Our  review  will  be  organized  as  follows; 

I.  Experimental  Evidence. 

A.  The  Solar  Corona 

B.  Solar  Corpuscular  Radiation  (SCR) 

C.  The  Heliomagnetic  Field  (HMF) 

D.  The  (MF-SCR  Interaction 

E.  Geomagnetic  Storms 

F.  Ring  Currents 

G.  The  Aurora  and  Airglow. 

II.  Solar  Corpuscvilar  Radiation 

A.  Kinetic  Properties  of  the  SCR 

B.  Parker's  Model 

C.  Chamberledn ' s  Model. 

III.  Geomagnetic  Storms 

A.  First  Riase  of  Storm 

1.  Early  Work 

2.  Cylindrical  Model 

3.  Plane  Model 

B.  Main  Phase  of  Storm 

1.  Ring  Currents 

2.  Aurorae. 

IV.  Steady-State  GMF-SCR  Intereu:tlon 

A.  One-Dimensional  Model 

B.  Two-Dimensional  Model 

C.  Three-DimensloneiL  Model 

D.  Analogies  from  Gasdyhainics 

E.  Stability. 
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I.  EXPERIMENTAL  EVIDENCE 


A.  THE  SOLAR  CORONA 

The  s\in  is  divided  into  three  regions:  (l)  the  photosphere,  which 
is  seen  with  the  naked  eye;  (2)  the  chromosphere,  visible  as  a  rosy  arc 
for  a  few  seconds  during  an  eclipse;  and  (3)  the  corona,  seen  as  a  pearly 
halo  surroiinding  the  lunar  disk  during  a  total  eclipse.  The  photosphere 
is  the  body  of  the  sun,  with  radius  R®  =  700,000  km  ;  the  corona  is 

the  solar  atmosphere,  and  the  chromosphere  forms  a  relatively  thin 
(~40,000  km)  transition  region  between  them  (AUer,  1953).  ^y  consider¬ 
ing  the  sun  to  be  a  blackbody  radiator  and  by  examining  the  energy  dis¬ 
tribution  of  its  radiation,  the  temperature  at  its  surface  is  found  to 
be  about  60OO  °K  , .  As  one  moves  radially  outward  through  the  chromosphere 

the  temperature  increases  to  1,000,000  -  2,000,000  °K  while' the  density 
16  0  3 

drops  from  10"^°  to  lo”  atoms/cm'^  at  the  base  of  the  corona. 

The  luminous  corona  is  highly  irregtilar  in  shape,  and  approximately 
one  to  two  solar  radii  in  width;  however,  rsuileLl  streamers  are  conspicuous, 
and  one  as  long  as  I5R®  (Chapman,  1961c)  has  been  recorded.  It  may  be 
more  useful  to  think  of  the  corona  as  continueLlly  expanding  outwards 
through  the  solar  system,  and  thus  constituting  the  SCR  as  well.  The 
luminous  corona  consists  of  three  components:  K,  E,  and  F:  the  K-com- 
ponent  has  a  continuous  spectrum  with  no  absorption  lines,  and  is  due  to 
scattering  by  free  electrcxis;  thus  its  llg^t  reaches  us  partly  polarized; 
and  from  measurements  of  polarization  and  luminosity  the  free  electron 
densities  are  calculated  (Van  de  Hulst,  1953J  Chapman  196lc). 

The  E-component  has  a  pure  emission  spectrum  of  bright  lines;  tem¬ 
peratures  determined  from  thermal  Doppler  widths  are  about  2  million  °K  ; 
whereas  the  degree  of  ionization  determined  from  the  emission  line  in¬ 
tensities  indicates  a  temperature  of  10^  ^  .  Other  determinations  also 

yield  temperatures  in  this  range.  The  differences  have  not  yet  'been  re¬ 
solved,  and  beyond  3R0  the  temperatures  can  only  be  inferred  indirectly 
from  electron  density  data  based  upon  eclipse  observations  (Chapman,  1961c). 
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Between  3R0  and  20R©  the  temperature  drops  from  10^  °K  to  3  x  10^  °K 
a  fact  which  leads  Parker  (see  II. B)  to  assume  the  corona  to  he  approxi¬ 
mately  idothermsG.  in  this  region. 

The  F-component  is  due  to  scattering  from  interplanetary  dust  between 
the  sun  and  the  earth. 

B.  THE  SOLAR  CORPUSCULAR  RADIATION 


Biemahn  (l95l).  demonstrated  that  radiation  pressure  alone  is 

ins\ifficient  to  account  for  the  observed  deflections  of  comet  tails 

passing  near  the  sun.  These  experienced  extreme  outward  accelerations 

which  he  believed  due  to  the  emission  of  a  neutrsJ.  stream  of  protons 

and  electrons  which  would  yield  densities  and  velocities  of  around 

100  protons/cm^  and  500  km/sec  at  1  a.u.  Kiepenheuer  (1953),  has 

sinmnarized  Biermann's  determinations  of  the  SCR  suid  the  work  has  been 

extended  by  Blermann  (l952>  1953>  1957#  19^0,  I961).  In  the  CEise  of  a 

If  #3 

solar  flare,  estimated  values  were  as  hig^  as  10  protons/ cm  and 
1500  km/sec. 

The  Explorer  X  rocket  probe,  launched  in  March,  196I,  measured 
quiescent  values  of  approximately  6-10  protons/cm^  and  300  km/sec 
(iCfy  Bull.  1962a;  Bridge,  et.  ed.,  1961).  During  the  experiment  a  solar 
flare  caused  a  geomagnetic  storm  at  the  eeurth;  this  appeared  to  Explorer  X 
as  a  svidden  rise  in  flux  density  and  peirticle  energies  beyond  37R  from 
the  geocenter.  Preliminary  results  Indicate  velocities  as  high  as 
700  km/sec  In  the  SCR  associated  with  the  flare. 

Explorer  X  measurements  further  Indicate  that  the  "temperature"  of 
the  SCR  Is  between  10^  and  10^  °K  (Bonettl  et  al,  1962;  Rossi,  I962). 
According  to  Bonettl,  et  al.  (I962)  the  lower  temperature  Is  suggested 
by  veurlatlons  In  particle  directions  ets  measured  by  Explorer  X;  the  upper 
teoqperature,  by  the  spread  In  particle  energies.  They  conclude:  '^These 
resvilts  are  rather  crude  but  they  Indicate  that  the  motion  of  the  par¬ 
ticles  Is  not  Isotropic.  Indeed  It  seems  that  the  thermal  motions  are 
predominantly  conflnedto  directions  parallel  to  the  magnetic  field  In 
the  frame  of  reference  moving  with  the  bulk  velocity. " 
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C.  THE  HELIOMAGNETIC  FIELD  (HMF) 


Until  the  Babcocks'  invention  of  the  solar  magnetograph  in  1952,  the 
HMF  of  the  undisturbed  sun  was  too  small  for  detection;  since  then  the  sun 
has  been  observed  to  have  a  dipole-like  field  of  approximately  1  gauss 
at  its  poles,  and  of  opposite  polarity  (Babcock  H.W.  and  Babcock,  H.D., 
1955)*  However,  towards  the  solar  equator  the  field  becomes  rather  dis¬ 
ordered,  although  still  approximately  one  gauss. 

The  streamers  in  the  corona  suggest  that  the  HMF  is  radial;  this  is 
confirmed  by  the  observations  of  Hewish  (1958)-  He  observed  the  scat¬ 
tering  of  radio  signails  from  the  Crab  nebula  as  it  moved  past  the  sun, 
and  found  the  HMF  to  be  radial  out  to  30R©  ,  which  was  the  limit  of 
his  observations.  Fiirthermore,  the  results  of  the  Pioneer  I  (Sonett, 
et.  al.,  1960a),  Pioneer  V  (Coleman  et.  sQ..,  I96O;  IGY  Bull.,  i960),  and 
Explorer  X  (IGY  Bull.,  1962a;  Bridge,  et,  al.,  1961)  space  probes  in¬ 
dicate  the  presence  of  an  interplanetary  field  of  approximately  2  -  IO7 
(17=  10  gauss)  .  Fields  of  the  order  of  6  -  127  were  measured  by 
Pioneer  I  and  Explorer  X  just  outside  the  magnetosphere,  and  Pioneer  V 
measured  fields  of  the  order  of  1-77  at  distances  of  five  million  km 
from  the  earth,  as  well  as  extreme  values  of  1^7  and  over,  which  were 
measured  simultaneously  with  the  onset  of  geomagnetic  storms  at  the  eeurth. 
Both  phenomena  were  correlated  with  the  appesirance  of  solar  flares. 

Pioneer  V  measurements  also  indicate  the  presence  of  a  steady  HMF  com¬ 
ponent  of  about  2.77  ,  a  figure  in  agreement  with  Parker's  theory  (see 

2 

II. B)  that  the  HMF  decreases  as  l/r  from  one  gauss  at  the  surface  of 
the  photosphere. 

The  study  of  rise  and  decay  times  of  cosmic-ray  activity  at  sea  level 
due  to  solar  flares  substantiates  the  notion  that  chcurged  particles  leaving 
the  sun  follow  curved  trajectories  along  magnetic  lines  of  force  originating 
in  the  sun,  and  do  not  come  directly  to  the  earth.  These  field  lines  are 
thought  to  stretch  outward  from  the  sun  with  a  characteristic  spiral  shape 
due  to  the  solar  rotation,  suggestive  of  water  Jets  from  a  rotating  lawn 
sprinkler  (Fig.  2).  Hence  the  name  "genrden  hose  effect."  The  persistence 
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HMF  and  effect  on  solar  cosmic  radiation. 


of  solar  cosmlc^^^  activity  long  erfter  the  flare  has  disappeared  indicates 
that  the  rays  may  he  trapped  within  the  HMF,  which  in  turn  is  stretched 
outward  due  to  the  pressure  of  the  trapped  particles. 

McCracken,.  (1962)  ,  has  described  the  cosmic-ray  flare  effects  and 
properties  of  the  HMF  which  may  be  deduced  from  them.  He  notes  that  the 
majority  of  events  are  due  to  solar  flares  on  the  west  limb  of  the  sion. 
Although  some  events  are  caused  by  flares  near  the  center  of  the  solar 
disk,  and  in  one  case  by  a  flare  in  the  eastern  of  the  disk,  rise 
times  of  solar  cosmic  ray  activity  measured  at  sea  level  (about  30  min¬ 
utes)  are  least  for  flares  occvirring  near  the  western  limb.  This  argues 
that  particles  diffuse  across  HMF  lines  as  they  travel  toward  the  earth, 
so  that  particles  from  the  western  limb  travel  the  most  direct  HMF  line 
from  sun  to  earth,  whereas  only  the  more  diffuse  edge  of  a  cloud  of 
particles  from  the  central  portion  of  the  sun  actually  reaches  the  earth. 
Furthermore,  McCracken  notes  that  very  energetic  particles  generally 
arrive  from  directions  making  angles  of  around  50°  with  the  earth-sun 
line;  this  agrees  with  the  direction  of  the  HMF  as  measxored  by  Explorer  X. 

The  Forbush  decrease  is  noted  Just  after  the  arrival  of  the  solar 
cosmic  rays:  this  is  a  decrease  in  the  amount  of  galactic  cosmic  radia¬ 
tion,  as  distinct  from  solar  cosmic  radiation,  which  increases.  Pioneer  V 
has  observed  a  Forbush  decrease  of  about  the  same  intensity  as  a  For¬ 
bush  decrease  measured  simultaneously  at  the  eeurth.  At  the  same  time 
the  average  HMF  increased  from  5  to  15  gammas  and  later  to  40  gammas, 
finally  returning  to  2  gammas  (aoiblent  value)  within  ^  hours  of  onset. 
This  lends  support  to  the  theory  that  the  HMF  is  conqpressed  and  blown 
outward  by  the  solar  flare  in  the  form  of  a  magnetohydrodynamic  shock 
wave,  and  that  the  increased  field  screens  out  the  normal  galactic  cosmic 
rays.  These  measurements  also  indicate  that  while  the  first  phase  of  a 
gecmagnetlc  storm  is  dependent  upon  solar  activity,  the  main  iflieise  is 
strictly  a  geoi^slcal  effect  whose  cause  lies  within  the  magnetosphere, 

^Cosmic  Rays  -  hl^ly  penetrating  corpuscular  radiation  originat¬ 
ing  outside  the  magnetosphere.  Solar  cosmic  rays  originate  in  the  stin; 
gedactlc  rays  originate  outside  of  the  solar  system.  Cosmic  Rays  con¬ 
sist  almost  entirely  of  positive  ions,  two- thirds  of  which  are  protons. 


since  the  main  phase  generally  persists  much  longer  than  48  hours. 

The  fact  that  the  solar  cosmic  radiation  continues  to  be  detected 
for  many  hours  and  even  days  after  onset  shows  that  the  Increased  fields 
also  keep  the  solar  radiation  from  moving  outward.  Parker  (1959)  l^as 
stiggested  that  this  may  also  be  due  to  the  presence  of  a  tangled  "shell" 
of  HMF  lines  which  always  enclose  the  inner  solar  system  somewhere  be¬ 
tween  Jupiter  and  Mars. 

D.  THE  GMP-SCR  INTERACTION 

The  Pioneer  I  space  probe  provided  the  first  direct  measurement  of 
the  distant  GMF  (Sonett,  et  al.,  1960a).  It  was  launched  on  October  11, 
1958,  one  of  the  quietest  days  on  record,  geomagnetically,  so  that  it 
may  provide  a  good  estimate  of  the  steady-state  fields,  as  well  as  a  lower 
bound  to  fluctuations.  Over  regions  from  3»7R  to  7R  >  and  frcmi  12. 3R 
to  14.8r  it  measured  the  component  of  total  GMF  perpendicular  to  its 
direction  of  motion;  however,  this  constituted  nearly  the  total  induc¬ 
tion  at  the  point.  The  authors  write  that  there  is  "close  genered  agree¬ 
ment  with  an  extrapolated  inverse-cube  law  within  10^  [fn  the  inner 
region] . "  In  the  outer  region  the  field  is  about  double  the  expected 
GMF,  euid  suddenly  drops  off  from  30  ganmas  to  around  6  gammas  at  I3.6R. 

The  relatively  sudden  drop  in  the  field  is  unmistakable,  but  the  resolution 
is  so  poor  that  one  could  only  say  that  the  transition  region  is  probably 
less  than  2000  km  thick,  i.e.,  less  than  2.51^  of  the  dimensions  of 
the  magnetos|diere.  The  direction  of  the  probe  was  close  to  local  noon 
while  traversing  the  outer  region. 

The  Explorer  X  probe  (IGT  Bull.  1962a),  measured  the  magnetic  field 
as  well  as  the  SCR.  At  6R  an  anguleu*  deviation  between  measured  and 
computed  fields  is  evident,  which  becomes  relatively  stable  at  a  dis¬ 
tance  of  llfl  ;  this  is  interpreted  as  indicating  an  appreciable  ex¬ 
traneous  field  superinposed  on  the  GMF.  This  effect  continues  out  to 
I9R  where  the  field  no  longer  shows  the  characteristic  l/r"^  decrease 
of  a  dipole  field.  The  authors  conclude  that  the  GMF  is  negligible 
here;  it  seems  more  likely  that  distortion  of  the  GMF  close  to  the 
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boundary  by  the  induced  field  have  altered  the  dipole  character  of  the 
GMF  beyond  recognition.  Between  20R  and  21. 5R  the  field  is  stable 
and  approximately  radial  from  the  sun  with  a  "stream  angle"  from  the 
radial  of  30°  to  40°  in  the  ecliptic  plane.  This  agrees  with  McCracken's 
figure  of  50°  (see  Fig.  2)  and  supports  the  notion  that  the  total  field 
at  the  interface  must  be  tangential  and,  therefore,  (at  this  position 
on  the  night  side)  nearly  parallel  to  the  SCR  and  the  HMF  embedded  in  it. 

The  actual  Interface  appears  to  be  at  21. 5R  ,  where  the  magnetic 
field  changed  dramatically  from  327  to  97  in  a  distance  of  less  than 
200  km.  There  was  also  a  violent  change  in  direction.  This  abrupt 
change  was  coincident  with  the  first  observations  of  SCR.  Thereafter, 
considerable  Irregxilar  fluctuations  as  wide  as  207  were  measured  out 
to  25R  after  which  they  settled  down  to  around  4-87  in  anqjlltude, 
superimposed  on  an  average  field  of  12  -  l6y  .  After  3^*5R  the  field 
declined  further  to  an  average  of  around  10 7  ,  interrupted  by  two  re¬ 
gions  in  which  the  fields  again  rose  to  around  20 7  ,  the  increase  being 
sustained  for  the  order  of  an  hour  or  more.  This  may  be  due  to  the  fact 
that  the  trajectory  is  nearly  parallel  to  the  expected  interface  in  this 
region  (Bonett,  et.  al.,  I962)  and  if  the  interface  is  slightly  unstable 
it  could  cross  and  re-cross  it  several  times,  and  apparently  does  re¬ 
enter  the  magnetosphere  finally,  Just  before  apogee  is  reached. 

The  information  obtained  should  eventually  yield  very  useful  inter¬ 
pretations  of  the  Interface  dimensions,  fields,  and  stability  of  the 
interface,  with  possible  implications  for  a  better  laiderstandlng  of  the 
Van  Allen  radiation  belts. 

E.  GEOMAGNETIC  STORMS 

Variations  in  the  GMF  sure  classified  as  S  (solar),  L  (lunar),  or 
D  (dlstvirbed),  with  numerous  subclassifications,  the  most  striking  of 
which  is  the  so-ceLLled  or  storm-time  variation.  This  is  the  under¬ 
lying  world-wide  systematic  variation  in  GMF  which  occurs  after  the  sud¬ 
den  ccsmnencement  (SC)  of  the  storm.  The  onset  of  the  storm  is  simultaneous 
to  within  one  or  two  minutes  edl  over  the  earth,  and  consists  mainly  of 
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a  sudden  sharp  rise  in  the  horizontal  (h)  component  of  the  GMP,  of  the 
order  of  20-3^7  >  which  may  persist  for  2  to  6  hours  and  is  known 
as  the  positive. initial,  or  first  phase  of  the  storm.  The  field  then 
drops  rapidly  below  its  original  value  euid  the  negative,  or  main  phase 
begins,  so-called  becuase  its  amplitude  is  generally  larger  than  that 
of  the  first  phase.  This  phase  may  last  for  several  days,  slowly  rising 
to  the  original  value  in  what  is  known  as  the  final,  or  recovery  phase 
(Mitra,  1952).  Changes  in  the  auroral  zones  are  generally  much  more 
severe  and  depend  strongly  upon  local  time. 

Verticeil  (Z)  component  changes  are  generally  much  smaller  than  the 
corresponding  H  changes,  and  of  opposite  sign.  Near  the  geomagnetic 
equator  Z  is  unaffected  by  storms;  however,  changes  in  H  are  greatest 
at  the  geomagnetic  equator. 

Storms  are  classified  as  moderate  i^Oy),  active  (15O7),  or  great 
(>2007)  ;  they  are  easily  distinguished  from  the  quiet  day  solar  varia¬ 
tion,  which  is  generally  of  the  order  of  +  IO7  .  The  veu'iatlon 

further  corroborates  the  concept  of  the  CMF-SCR  interaction:  it  is 
greatest  during  the  day  and  in  the  summer  hemisphere;  it  is  50  -  100^ 
greater  at  sunspot  maximum  than  at  sunspot  minimum  (Mitra,  1952). 

Storms  have  been  definitely  correlated  with  the  appearance  of  soleir 
flares,  although  there  are  numerous  cases  in  which  fairly  strong  storms 
have  occurred  without  any  visible  solar  activity.  Storms’ are  also 
acccjnpanied  by  ionospheric  disturbances  and  polar  auroral  activity. 

F.  RING  CURRENTS 

The  existence  of  a  ring  current  high  above  the  earth  was  postulated 
by  Chapman  and  Ferraro  (1933)  "to  explain  the  main  phase  of  the  geomagnetic 
storm,  which  could  not  be  accounted  for  by  SCR  compressing  the  GMP,  as 
in  the  first  phase.  Ry  means  of  a  harmonic  analysis  of  field  fluctuations 
in  storms  they  showed  that  the  main  phase  may  be  ascribed  to  a  west- 
ward-flowing  ring  current  circling  the  esurth  at  a  height  of  a  few  hundred 
kllcmieters,  with  intensities  of  the  order  of  100,000  an^res,  depending 
upon  the  strength  of  the  storm  (Mitra,  1952 ) .  Later  estimates  placed  the 
altitude  of  the  current  at  around  4r  to  7R  . 
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The  first  direct  measurements  of  this  current  were  those  of  Explorer  VI 
and  Pioneer  V  rocket  probes  (Smith,  et.  aJ..,  I96O;  Sonett,  et.al.,  1960b), 
which  discovered  anomalies  in  the  GMP  between  4r  and  8R  ,  explainable 
on  the  assumption  of  a  ring  current.  With  its  geometry  and  intensity 
determined  to  give  closest  agreement  with  the  data,  the  authors  found 
a  toroidal,  westward-flowing  ring  current  of  radius  9R,  circular  cross- 
section  of  radius  3R#  and  current  of  the  order  of  5  Megamps.  The  eixis 
of  the  ring  was  taken  parailel  to  the  earth's  axis. 

It  should  be  noted,  however,  that  the  model  was  constructed  on  the 
basis  of  very  little,  and  inconclusive,  evidence  which  may  only  be  con¬ 
sidered  strongly  sviggestlve,  at  best.  Later  probes.  Explorers  X  and  XII, 
foiond  no  evidence  of  ring  cvirrent.  Nonetheless,  it  appears  that  the 
concept  of  a  ring  ctirrent  giviiag  rise  to  the  main  phase  of  the  storm 
has  much  to  commend  it,  although  its  existence  and  source  remain  in 
doubt . 

G.  THE  AURORA  AND  AIR  GLOW 

We  shall  here  review  only  those  auroral  phenomena  which  seem  to  be  most 
directly  related  to  the  OlF-SCR  interaction.  The  most  important  of  these 
is  the  well-established  correlation  between  auroral,  magnetic,  and  solar 
activity,  the  eleven-year  cycle  being  cleeurly  evident  in  all  three 
I^enomena. 

Auroras  generally  occur  in  hl^  latltixdes;  contours  of  equal  fre¬ 
quency  of  occurrence,  or  isochasms,  are  rou^ly  concentric  ovals,  with 
centers  approximately  at  the  magnetic  pole.  The  maximum- frequency 
Isochasm  in  the  northern  hemisphere  lies  along  67°  north  magnetic  latitude; 
a  similar,  but  less  well-known  phenomenon  occurs  In  the  southern  hemi¬ 
sphere.  At  this  latitude  the  awora  averages  243  days/year,  decreasing 
to  around  200/yr  at  72%  and  lOO/yr  at  62%  magnetic  latitude  (Chapman, 

1961a). 

Riotoelectrlc  observations  of  aurorel  spectra  during  the  day  show 
that  activity  Is  largely  restricted  to  local  midnight  +  3  hours.  Using 
the  dipole  as  a  reference  axis  It  is  found  that  the  prlnclpeiL  maximum 
of  auroral  Intensity  occurs  about  one  hour  before  magnetic  midnight. 
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Althou^  the  aurored.  spectrm  is  made  up  largely  of  oxygen  and  nitrogen 
lines  exhibiting  low-temperature  thermal  Doppler  broadening,  a  smell  part 
of  the  li^t  comes  from  hydrogen  QJ  and  P  lines .  When  viewed  sideways 
these  lines  show  Doppler  broadening,  indicating  thermel  velocities  of 
the  order  of  400  km/sec  ;  when  viewed  from  below,  along  a  GMF  line 
which  the  aurora  tends  to  follow,  there  is  an  electromagnetic  Doppler 
shift  of  the  entire  spectrum  indicating  that  the  emitting  particles  are 
descending  with  mean  speeds  of  around  ^00  km/sec  ,  and  some  as  high 
as  3000  km/sec  (Chapman,  196la).  This  is  comparable  to  the  SCR  drift 
velocity  of  300  km/sec. 

We  have  already  mentioned  the  correlation  between  auroras  and 
geomagnetic  stormsj  in  IGY  Bull.  No.  45,  (l96la)  we  find  an  instance  of 
direct  observation  of  the  "dximping"  of  particles  from  the  inner  Van  Allen 
belt  into  the  auroral  zones  diirlng  a  geomagnetic  storm. 

The  night  "airglow"  or  "ni^tglow”  acco\ints  for  50^t  of  the  lum¬ 
inescence  of  the  night  sky,  apart  from  auroras.  It  is  a  sort  of  per¬ 
manent  aurora,  and  always  present;  althou^  its  intensity  is  small  com¬ 
pared  to  auroral  Intensities,  the  question  of  the  precise  distinction 
between  aurora  and  airglow  is  not  yet  resolved  (Chamberlain,  1961b). 

Thus  far  the  airglow  has  only  been  observed  at  night;  Intensities 
eure  hl^er  at  twilight  and  before  dawn  than  at  mldnlf^t  minimum  by  a 
factor  of  around  fifteen.  The  "dayglow”  spectrum  is  believed  similar 
to  that  of  the  "twillghtglow”  with  intensity  as  great  or  greater;  however, 
attempted  rocket  observations  of  dayglow  have  resulted  in  contradictory 
estimates  (Chamberlain,  1961b).  In  addition  to  this  quasi-regular  diurnal 
behavior,  there  is  also  a  seasonal  variation,  and  Irregular  behavior 
which  is  closely  correlated  with  magnetic  activity. 


II.  SOLAR  CORPUSCULAR  RADIATION 


A.  KINETIC  PROPERTIES  OF  THE  SCR 

The  following  calctilations  of  velocity,  mean  free  path,  and  conduc-^ 
tivity,  are  "based  on  Delcroix  (i960)  Chap.  11,  in  which  it  is  shown  that 
for  a  strongly  ionized  gas  an  effective  collision  frequency  is  defined  as 

V  =  0.0876  nT’^/^log^A  ,  (1) 

in  the  case  of  proton-proton  collisions.  The  temperature  of  the  tons  is 
T,  n  is  their  numher/cm^,  andA  is  the  average  ratio  of  Dehye  length  to 
distance  of  closest  approach  for  proton-proton  Coulomb  collisions.  If 
we  take  the  Explorer  X  figures  for  the  SCR  at  1  a.u.  (astronomical  unit) 
of  n  =  lO/cm^  ,  T  =  500,000  °K  ,  then 

log^A  =  27.9  j  V  *  6.90  X  10"®  sec”^  ,  (2) 

and  the  electrical  conductivity  is  (Spitzer,  1956,  Chap.  5) 

c  =  1.9^  X  10^  mhos/meter  ,  (3) 

considering  ion-electron  and  electron-electron  collisions  in  a  singly-ionized 
plasma.  The  mean  thermal  velocity  of  eua  ion  is 

v^  =  (Skr/rtm)^^^  =  103.4  km/sec,,  (4) 

and  the  mean  free  path  is 

«  v^/v  *=  15.6  X  10®  km  =  2.35  X  10^  R  *  10.0  a.u.  (5) 
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Even  with  a  more  conservative  estimate  of  T 


=  1C)^.°K  we;  find 


•  V  a  7«0T  x  lO'T  sec  ^  v^  =,  46.3  km/sec 

s  6.54  X  km  =  0.44  a.u.  ,  (6) 

so  that  the  mean  free  path  of  SCR  ions  at  the  earth's  orbit  is  considerably 

larger  than  the  geomagnetic  cavity  covild  possibly  be,  even  larger  than 

the  dimensions  of  the  inner  solar  system. 

One  may  apply  the  same  formula  at  the  base  of  the  solar  corona,  where, 
6  n  ft 

let  \is  say,  T  =  10  a  ,  n  =  10  /cm"^  ;  then 

V  =  0.18  sec  v^  =  l46.4  kn/sec;  =  8l4  km;  ;  (T) 

this  last  figure  may  be  con^iared  with  the  solar  scale  hei^t 

h  =  kT/mgg  =  31,000  km  ,  (8) 


where  g  is  the  solar  gravity. 

8 

We  find  the  thermal  pressure  of  the  SCR  at  1  a.u.  is 


for 


P 


2nkr 


10^  °K  , 

■  (4 


X  lO”^^  dynes/cm^  , 


(8.1) 

(8.2) 


and 


2  nmV^  «  3  X  10  ®  dynes/cm^  , 


(8.3) 


the  pressure  of  the  SCR  reflected  from  the  front  of  geomagnetic  cavity. 
Thus,  on  the  dayslde  the  thermal  pressures  are  generally  neglected  com¬ 
pared  to  the  dynamic  pressure,  2nmV^.  . 
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B.  PARKER'S  MODEL 


Inspired  by  Chapman's  (l95T#  1959)  hydrostatic-conductive-equllibrlvim 
model' of  the  corona,  Parker  (l958a#l3#  1959^  1960a,b,  196la,b,c)  has 
constructed  a  hydrodynamic  continuum  model  of  the  expanding  corona  with 
which  he  has  attempted  to  deduce  observed  coronal  temperatures,  densities, 
and  expansion  velocities.  The  best  review  of  his  work  is  found  in  Parker. 
(1961a), although  it  omits  many  of  the  mathematical  detaila. 

The  basic  ingredients  of  the  continuum,  solutions  attempted  by  Chapmeui, 
Parker,  and  later  Chamberlain  are: 

(1)  neglect  electromagnetic  fields  in  the  corona  apd  consider  it 
to  be  a  quasi-neutral  single  fluid; 

(2)  assume  spherical  symmetry;  neglect  solar  rotation; 

(3)  conservation  of  matter  gives  a  simple  integrated  continuity 
relation  in  sphericeLl  symmetry, 

2  2 

nr  w  =  n^a  w^  =  const.  ,  (9) 

where  n  =  number  density  of  ions,  w  is  the  drift  velocity  of  the 
expanding  corona,  r  is  the  heliocentric  distance  from  the  center 
of  the  sun,  and  subscripts  (q)  refer  to  values  at  r  =  a,  the 
"base"  of  the  corona; 

(4)  assume  an  isotropic  pressure  tensor,  although  it  may  be  ex¬ 
pected  that  the  high  velocity  of  the  stream  and  the  presence  of 
the  HMF  will  eventually  introduce  anisotropies: 

(^)  assume  an  equation  of  state  for  the  corona, 

p  =  2  nkT  (perfect  gas)  ; 

(6)  assume  validity  of  the  classical  single-fluid  Navler-Stokes 
equation 
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S  w/St  +  (w:grad)v  =  -  (l/nm)grad  P  -  §3  ,  ; 

where  p  =  pressure,  m  =  ion  mass,  is  acceleration  due  to 

solar  gravity; 

(7)  assume  some  form  of  the  heat  conduction  equation  to  apply t 

div  (K  grad  T)  =  Bo/St  , 

where  K  is  the  thermal  conductivity  and  Q  is  the  energy  density; 

(8)  consider  only  the  steady-state,  h/ht  =  0. 

fTViapmAr)  neglected  the  convective  term  in  (6),  and  Parker  noted  that 
this  results  in  a  pressure  which  approaches  '  10  dy/cm  at  r  =  “  , 

much  larger  than  the  interstellar  pressure  of  10  dy/cm  thought  to 
actually  exist  there.  Thus,  the  corona  cannot  he  in  static  and  con¬ 
ductive  equilibrium,  Parker's  approach  includes  the  convective  term 
in  (6),  but  avoids  (7)  by  assuming,  a  priori,  that  for  some  undetermined 
r  =  b  ,  there  is  a  region  a  <  r  <  b  over  which 'the  corona  is  iso¬ 
thermal,  while  its  temperatxu'e  decreases  adiabatically  for  r  >  b  . 

Thus  (7)  and  (8)  are  not  satisfied,  and  there  is  no  conductive  equilibrium. 
For  a  =  10^^  cm  {1.^3  R©)  Parker  sets  T^  =  10^  "^0  “  ^  ^  loVcm^ 

In  his  analysis  Parker  derives  a  family  of  solutions,  all  of  which 
are  either  physically  impossible,  or  else  yield  very  low  drift  velocities 
of  the  order  of  30  km/sec  ,  hardly  the  "solar  wind"  predicted  by 
Biermann,  but  more  like  Chamberlain's  ’solar  breeze.  However,  it  is 
Parker's  contention  that  there  is  one  and  only  one  "singular"  solution 
which  behaves  properly;  i.e.,  gives  a  physically  meaningful  solution 
with  hi^  velocity  far  from  the  sun  and  low  velocity  near  the  sun.  This 
is  the  crucial  issue,  because  Chamberl8d.n  (i960,  196la)  claims  that 
Parker  has  Integrated  his  equations  Incorrectly,  that  his  slng\ilar 
solution  is  not  mathematically  valid.  Chamberlain  believes  that  a 
solution  is  either  high-  or  low-velocity,  and  cannot  'cross  over;  that 
Parker's  solution  implies  finite  energy  per  particle  at  infinity,  which 
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in  t\irn  must  lii5)ly  mass  transport  of  energy  so  large  as  to  make  heat 
transfer  by  conduction  negligible,  so  that  the  adiabatic  approximation 
should  hold  almost  to  the  corona,  contradicting  Parker's  assumption  of 
an  isothermal  corona  between  a  and  b.  Furthermore,  to  maintain  coronal 
temperature  out  to  1OR0  or  more  requires  too  much  energy.  Despite  re¬ 
buttal  and  counter-rebuttal  the  controversy  over  these  points  remains 
in  force  at  the  present  time  (1962). 

Althoiagh  b  is  at  Parker's  disposal,  is  not,  due  to  the  highly 

"singular"  nature  of  his  solution,  and  the  assimied  value  of  -.'T^  Cal¬ 

culations  of  n,  V  for  several  combinations  of  T  and  b  (Parker,  1960b, 

^  u 

pp.  856-858),  yield  T  *  10  °K.  ,  b  =  8a  =  11.4  R©  ,  and 
n  =  20/ cm^  ,  w  =  310  km/sec  at  1  a.u.  ,  which  is  in  very  close 
agreement  with  the  Explorer  X  space  probe  results  (see  I.D.)  and  was 
published  six  months  before  the  laiinching  of  that  probe.  In  addition, 
it  should  be  noted  that  for  =  1.22  x  10^  °K  ,  Parker's  model 

of  the  isothermal  corona  exhibits  agreement  with  density  observations 
based  on  electron  scattering  which  is,  in  his  words,  "...so  extraordinary 
as  certainly  to  be  fortuitous."  But  it  is  very  suggestive,  nonetheless. 

However,  there  are  numerous  difficulties,  including  the  fact  that 
by  using  Parker's  figwes,  Tq  =  1.22  x  10^  °K  ,  and  n^  =  2  x  10^/cm^ 
his  adiabatic  temperatvire  relation, 

(T/Tq)  *  (n/n^)^/3  >  (10) 

yields  T  =  16,000  °K  for  n  =  1  a-u.  A  figure  which ,  . 

may  be  coiQ>ared  with  Chamberlain's  (15-20,000  °K)  ,  Chapman's 

(440,000  °K)  and  the  Explorer  X  results  (lO^  -  10^  °K)  .  The  similarity 

to  Chamberledn,  who  also  uses  a  hydrodynamic  approach,  and  the  disagreement 
with  empirical  estimates  coupled  with  the  electron-density  agreement  near 
the  sun,  and  with  the  values  of  ion  mean  free  path  calculated  in  the  pre¬ 
vious  section  (II.A),  Indicate  that  the  cleisslcal  continuum'  approach  may 
Indeed  be  veLLld  near  the  sun,  but  that  the  large  mean  free  paths  at  1  a.u. 
require  that  we  go  over  to  a  free-molecule  flow  at  some  distance  between 
the  sun  and  the  earth. 
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Parker  derived  an  expression  for  the  HMF  by  assuming  that  particles 
leave  the  sun  with  constant  velocity  and  drag  the  HMF  lines  with  them. 

Thus,  in  the  system  rotating  with  the  sun  these  particles  describe  spiral 
arcs  (the  "garden-hose"  effect)  and  the  HMF  must  be  tangent  to  these  arcs. 

By  combining  this  requirement  with  div  B  =  Q  ’  ,  h’e  finds  field  linds  in 
the  plane  of  the  ecliptic  being  stretched  out  and  bent  by  the  SCR,  with 
radial  component  depending  on  l/r^  ,  the  azimuthal  component  on  l/r  . 
Parker.  (1959),  states:  "The  field  at  the  orbit  of  the  earth  is  very  nesirly 
radial....  The  field  is  unstable  because  of  the  anisotropic  expansion  of 
the  gas  as  it  canes  out  from  the  sun.  The  collision  rate  is  extremely 
lowj  hence  if  the  gas  expands  in  directions  perpendicular  to  the  radius, 
but  not  in  the  radial  direction,  its  thermal  motions  become  anisotropic 
and  the  gas  Itself  becomes  unstable.  The  field  becomes  disordered."  This, 
in  Parker's  view,  forms  the  magnetic  shell  which  he  feels  is  responsible 
for  the  diffusive  decay  of  cosmic  rays.  The  Forbush  decrease  in  gEJLactic 
cosmic  rays  during  a  geomagnetic  storm  is  explained  by  the  fact  that  storm 
particles  are  more  energetic,  their  stream  lines  have  less  curvature  and 
stretch  out  the  HMF  forming  a  magnetic  shock  wave  (McCracken,  I962). 

C.  CHAMBERLm'S  MODEL 

Initially,  Chamberlain,  (i960),  applied  Jeans'  evaporative  theory  to 
the  calculation  of  the  expanding  solar  corona,  hud  more  recently  (Chamber- 
lain,  1961a),  he  has  adopted  Parker's  approach,  but  with  a  fundamental 
difference:  Chamberlain  has  used  a  heat  equation  derived  from  the  first 
law  of  thermodynamics 


dq  =  du  +  pdv  ,  (U) 

where  q,  u,  v  eure  heat  absorbed.  Internal  energy,  and  volume,  respectively, 
per  unit  mass.  Since  the  heat  absorbed  per  unit  volume  per  second  Is 
given  by  Fourier's  heat-conduction  equation  to  be  div  (K  grad  T)  ,  then 
the  heat  absorbed  per  gram-sepond  is  (Llepmann  suid  Roshko,  1957^  Chap  I3): 
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(l/mn)  div  (Kgrad  T)  =  dq/dt  =  du/dt  +  p  dv/dt  .  (12) 


If  the  entire  system  is  moving  we  may  replace  the  d/dt  operator  with  the 
mobile  operator  S/St  +  w*  grad'',  measured  hy  the  stationary  observer. 

Since  u  =  31^/2m  and  v  =  l/mn  ,  we  set  S/St  =  0  and  find: 

div  (K  grad  T)  =  (3nk/2)  w*  grad  T  -  liN.  grad  n  ,  (13) 

where  the  gas  is  understood  to  be  monatomic.  Thus,  if  he  taJces  K  =  ^0^^^ 
Chamberlain  finds  representative  vedues  of  <  21,500  °K  and  n  =  30/cc 
at  one  avu.  .■  These  are  comparable  to  the  values  obtained  by  the  evapora¬ 
tive  model,  but  far  from  the  empirical  estimate  of  10^  -  10^  °K  .  How¬ 
ever,  Chamberlain  has  initially  assumed  the  vanishing  of  particle  expansion 
energies  at  <»>,  so  that  his  "breeze”  conclusion  is  to  be  expected,  and,  of 
itself,  does  not  preclude  a  solar  wind,  only  its  inevitability  as  the  re¬ 
sult  of  Parker's  "singular"  solution.  This  paper  is  to  be  reccsnmended  for 
its  clarity. 

Although  the  ion  mean  free  path  at  1  a.u.  is  of  the  order  of  an  astron¬ 
omical  unit,  it  may  not  necessarily  invalidate  a  continuum  approach,  since 
the  lArmore  radius  of  the  protons  in  the  SCR  is  3000  km  at.  l.a.u,  for 
a  HMP  of  I7  ,  a  relatively  small  length.  Thus  peurticles  revolve  about 
field  lines  which  transmit  the  stresses  of  the  coronal  expansion  throu^- 
out  the  inner  solar  system.  Thus,  one  may  perhaps  treat  the  problem  near 
1  a.u.  ,  by  means  of  an  anisotropic  pressure  tensor,  as  suggested  by 
Bonetti,  et.  al.,  (1962.) 
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III.  GEOMAGNETIC  STORMS 


A.  FIRST  PHASE 

1.  Early  Work 

At  first  glance  it  may  appear  out  of  order  to  discuss  the  nonsteady 
GMF-SCR  interaction  before  investigating  the  steady-state;  hut,  histori¬ 
cally,  it  was  the  nonsteady  phenomena  which  first  attracted  the  attention 
of  science,  and  eventually  led,  within  the  last  decade,  to  the  concept  of 
the  steady-state  interaction.  In  I892,  Lord  Kelvin  demonstrated  that  the 
storms  could  not  he  directly  due  to  variations  in  the  HMF,  and  Hckle's 
measurements  confirmed  this  fact  so  that  the  interaction  had  to  he  either 
corpuscular  or  radiative  (ultra-violet).  In  I896,  Blrkeland  showed  that 
electron  beams  shot  out  of  the  svua  would  he  deflected  by  the-  GMF,  by  pro¬ 
jecting  cathode  rays  toward  a  magnetized  sphere  with  the  exotic  name  of 
Terrella.  Ingpact  was  generally  in  the  polar  regions,  (Mltra,  1952,  p.  464), 
and  since  then  the  theories  of  geomagnetic  storms  and  auroral  i^encnena 
have  been  linked  together. 

Following  this,  Stormer  developed  the  theory  of  orbits  of  Individual 
charged  particles  In  a  magnetic  dipole  field.  Ry  1931#  Chapman  eud  Ferraro 
(193181^  P*  80)  were  able  to  state  that  "...our  work  confirmed  Llndemann’s 
conclusion  that  the  only  admissible  kind  of  stream  Is  one  that  la  electro¬ 
statically  neutral  to  a  very  high  degree  of  approximation.”  Mltra  (1952, 
Chap.  IX)  offers  an  an^le  review  of  the  work  of  Blrkeland  and  Stormer) 
however,  the  emswers  to  the  various  problems  discussed  In  this  review  are 
not  to  be  found  In  the  realm  of  classical  slngle-pcurtlcle  physics,  although 
such  methods  hove  afforded  us  numerous  valuable  Insights. 

Most  of  the  theoretical  foundation  for  present  thought  qn  the  sub¬ 
ject  of  the  GMF-SCR  Interaction  Is  based  on  the  work  of  Chapman  and  Ferraro, 
collectively  and  Individually,  In  which  they  attempted  to  expledn  the  origin 
and  behavior  of  magnetic  storms  by  means  of  certain  mathematical  IdeeiLlza- 
tlons.  The  first  such  was  the  application  of  Maxwell's  solution  (Maxwell, 
1881),  for  a  magnetic  dipole  moving  parallel  to  a  perfectly  conducting 
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plane  sheet.  This  showed  the  lines  of  force  to  he  parallel  to  the  boundary 
except  at  two. neutral  points,  N,  where  the  field  vanishes  (see  Fig.  3). 

The  image  dipole  doubles  the  field  component  parallel  to  the  boundary  and 
concels  the  normal  component,  but  it  vanishes  within  the  conductor. 

Chapman  and  Ferraro  tried  to  demonstrate  the  existence  of  the  geomag¬ 
netic  cavity  and  ascertain  its  size  along  the  earth- sun  line  by  consider¬ 
ing  the  forces  on  a  "surface  layer"  due  to  accumulation  of  matter,  change 
of  momentum  of  matter  entering  the  layer,  and  magnetic  stresses  fran  the 
axogmented  GMF  inside  the  cavity.  The  results  obtained  then,  by  laborious 
approximation,  agree  very  closely  with  more  recent  results  obtained  by 
Ferraro  (1952,  Section  5.6)  with  the  aid  of  an  electronic  computer. 
However,  both  treatments  of  this  difficult,  time-dependent  problem  fail 
to  account  satisfactorily  for  continuity  relations  at  the  surface  (see 
IV. A.)  and  give  rise  to  a  very  fundaments^,  contradiction  concerning  mo¬ 
mentum  balance  across  the  surface  layer.  However,  their  concept  of  the 
geomagnetic  cavity  has  been  borne  out  by  recent  measurements. 

As  a  result  of  their  investigations  Chapmeui  and  Ferrsuro  developed  a 
model  of  the  GMF-SCR  interaction  which  also  embodied  the  germs  of  the  ring 
current  idea,  and  a  possible  explanation  of  the  correlation  with  aurorsiL 
activity.  They  conclxided,  (1931a,  p.83):  "in  the  csise  of  a  stream  from 
the  Sun  advancing  towards  the  Eeud^h,  a  hollow  space  round  the  Eeurth  is 
formed  in  the  stream. . .open  at  the  back  of  the  Eeurth. . .part  of  the  stream 
which  has  collected  near  the  Esurth  rems^ns  for  a  time,  probably  in  the 
form  of  a  ring  around  the  equator;  this  gradusHy  dlsappesurs  by  the 
passage  of  the  ions  and  electrons  along  the  Earth's  lines  of  force,  into 
the  atmosphere  in  high  latitudes." 

2.  The  Cylindrical  Model 

.  In,  19^.  Chaimah  and  Ferraro  preaeated  the  cyllndei*  mpdeily  .Ih,  wbl^  a 
cylindrical  sheet  of  plqdsma  collapses  from  infinity  upon  a  coaxial  mag¬ 
netic  field  (see  Fig.  4).  The  original  (permanent)  field  is 


FIG.  4— The  Cylinder  Model. 
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(n>2)  Hp  =  H^CB/rA  Ap  =  (-  H^rAii  -  2))  (B/r)*^"^ 


W 


vhere  H  is  in  the  z-direction,  and  its  vector  potential,  A  ,  is  in  the 

p  .  '  P 

♦■^direction,  tangential,  to  the  cylinder.  The  sources  of  are  taken 
to  be  COT^letely  unaffected  by  the  passage  of  the  current  sheet. 

In  addition  to  Maxwell's  eq.vtations,  the  authors  employ  the  equations 
of  motion  for  the  individual  particles  subject  to  the  Lorentz  force.  The 
fields  are  derived  from  an  electric  scEG.ar  potentisd.  ♦  =  ^(r)  and  a 

magnetic  vector  potential  A  =  A(r)  in  the  0-direction.  The  ♦  is  to 
account  for  electrostatic  forces  due  to  possible  charge  separation;  the 
vector  potential  A  incl\ides  the  effect  of  ion  and  electron  currents  in  the 
cylinder. 

Althovi^  the  authors'  calculations  are  straightforward,  they  are 
involved;  a  simpler  calculation  which  makes  use  of  their  insights  suffices 
to  reproduce  the  mein  features  of  their  aneilysis.  The  electromagnetic 
Hamiltonian  of  a  particle  of  charge  q,  mass  m,  is: 


»  (l/2m)  (£  -  q^A/c)^  +  q^  =  (m/2)  (u^  +  v^)  +  q^  =  mU^/2  ,  (ip) 

where  the  peurticle  velocity  In  (r,  0,  z)  -  coordinates  is  v  =  (u,  v,0)  , 
and  at  infinity,  v  *  (U,  0,  O)  .  Since  ♦  =  ♦(r)  ,  the  0-coniponent 

of  the  canonical  momentum  is  a  constant,  zero  in  this  case  because  A  and 
V  vanish  faster  than  l/r  at  infinity;  therefore, 

mv  =  -  qA/c  ,  (l6) 


.  and  thus, 


m  V 
e  e 


electrons 

ions 


(17) 


l.e.,  the  anguleur  momentum  of  the  cylinder  is  conserved,  which  might  have 
been  anticipated  on  grounds  of  symmetry. 
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On  the  basis  of  Eq,.  (17)#  one  may  assume  the  current  J  to  be  due  to 
the  electrons, 


J  =  Jg  =  -Qv/2ncr  ;  =  -  Qy^/c  ,  (18) 


where  Q  represents  the  amount  of  positive  charge  per  unit  length,  and 
is  the  magnetic  vector  potential  due  to  the  electrons  evaluated  just 
inside  the  surface  of  the  cylinder. 

From  Eq,.  (16),  we  have 


m  V 
e  e 


=  (e/c)  (Ap  +  A^) 


(19) 


Substituting  for  A^  from  Eq.  (18)  yields 


m^v^  =  -  [cm^R^Ho/Cn  -  2)  aQr"  '  ,  (20) 


where  a  =  1  +  m^c  /eQ  ;  (Jhapman  (i960)  shows  that  a  =  1  for  eOl 
cases  of  geomagnetic  interest.  The  Induced  field  inside  the  cylinder  is 

He  ?  4«Jg  »  [2H p/a(n  -  2)]  .  (2l) 

One  may  now  consider  the  Inward  radial  motion  of  the  cylinder  as 
a  whole,  the  forces  upon  it  given  by  the  difference  in  magnetic  pressure 
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across  the  surface.  By  symmetry,  u  =  u(t)  =  u[r(t)]  ,  so  that  the 
equation  of  motion  per  unit  length  of  cylinder  is  (a  =  1  ,  n  =  3)  ‘ 

(Qm  /e)  du/dt  =  (Qm7e)u  du/dr  =  2jtr(H|  +  2H  H  )/8*  .  (22l) 

i  ®  P 

Since  H  =  2H  ,  integration  yields 
6  P 


-  u^  =  (el^  (R/r)^  ,  (23) 

and  setting  u  =  0,  we  find  the  stopping  radius,  r^  : 

(r^/Rf  =  (e/Qmj^)^/^  (HqR/u)  .  (24) 

The  total  Increase  in  magnetic  energy  per  unit  length  at  the  stopping 
radius,  R  =  r^  ,  is 

^\ag  '  (’«-0/8«)  C(HpQ  +  H^q)^  -  “  (Q^iU^/e)  ,  (25) 

or  twice  the  available  kinetic  enei^.  Chapman,  (i960,  p.  924)  .  has 
noted  this  "apparent  peuradox"  which  he  says  is  due  to  postulating  that 
the  permanent  field  is  unaffected  by  the  passage  of  the  cylinder,  so  their 
mutviEG.  Inductance  is  zero,  and  cxie  considers  only  the  Induced  field; 
i.e.. 


-  (■•B/8)  4o  -  ('»"7/2e)  • 

However,  it  is  probably  more  correct  to  say  that  ■•■he  additional  energy 
is  supplied  by  the  sources  of  the  field  in  order  to  maintain  the  perma¬ 
nent  field  during  the  passage  of  the  cylinder. 
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It  might  be  possible  to  test  this  theory  in  its  essentials  by 
covering  a  thin  metal  cylinder  vith  a  very  uniform  shaped  charge  of 
explosive,  and  placing  a  uniform  magnetic  field  inside  by  means  of  a 
large  magnet  and  a  search  coil  at  its  center  to  measure  5B/5t  .  The 

explosion  should  vaporize  the  metal  and  drive  it  inward,  giving  rise 
to  a  situation  much  like  that  treated  by  CHiapman  and  Ferraro. 

As  an  approach  to  the  real  physical  situation,  the  shortcomings  of 
this  model  are  evident;  however,  it  constituted  an  important  step  for¬ 
ward  in  establishing  the  current  view  that  the  SCR  is  retarded  and  re¬ 
flected  from  the  (MF  as  well  as  providing  some  theoretical  background 
for  Ferraro's  1952  paper.  It  edso  provided  the  clear  indication  that 
the  size  of  the  geomagnetic  cavity  would  be  determined  by  a  pressure 
balance  across  its  boundary,  involving  the  magnetic  field  on  one  side 
and  the  change  in  peurticle  momentum  on  the  other.  The  concept  of  "specular" 
reflection  of  particles  arises  simply  as  a  consequence  of  the  symmetry: 
since  the  potentials  involve  only  one  coordinate,  r,  the  trajectories 
of  indlvidueLl  particles  must  be  symmetrical  with  respect  to  their  perigee 
radius,  giving  the  effect,  at  large  distances,  of  specular  reflection. 

More  recently,  Chapnan  and  Kendall  (1961b);  have  made  a  calculation 
of  the  motion  of  the  collapsing  cylinder  in  the  case  that  the  cylinder 
particles  8u:e  projected  obliquely  inward;  rather  than  radially.  The 
purpose  was  to  obtain  some  idea  of  the  shape  of  the  GMF-SCR  Interface 
in  the  equatorial  plane. 

3.  The  Plane  Model 


In  19^2,  Ferraro  published  a  plane  model  of  a  geomagnetic  storm. in 
which  an  infinite  stream  of  plasma  at  T  *  0^  advances  in  the  (-z)' 
direction  into  a  permanent  magnetic  field  •  Hq(r/z)^  directed 

In  the  y-dlrectlon  (see  Fig;  5)«  ^  assuming  a  total  Vector  potential 

A  «  A(z)  In  the  x-dlrectlon,  and  an  electrostatic  field  derivable 
from  a  scalar  potential  t  «  ^(z,t)  due  to  charge  sepetratlon,  Ferraro 
derives  approximate  solutions  from  the  time-dependent  slngle-peurtlcle 
Lagranglans  In  a  manner  similar,  In  many  respects,  to  the  solution  of 
the  Cylindrical  Model.  However,  since  the  stream  Is  continuous,  and 
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infinite,  and  the  stream  front  is  in  motion,  the  time -dependence  greatly 
complicates  matters. 

Ferraro  derives  the  equation 


=  43t(e/c)^  (n^/m^  +  n^/m^)A  ,  (26) 


and  after  ass\jming  n^  =  n^  =  n  he  derives  an  approximate  repre¬ 


sentation  for  A  hy  the  WKB  method: 


A  =  0!  X  exp  [-  /X  dz]  ;  =  4rtne^/m  c^ 

^0 


(27) 


3y  proceeding  from  this,  he  is  then  able  to  demonstrate  that  X  =  const., 
vith  a  possible  change  of  the  order  of  unity  right  at  the  surface.  How¬ 
ever  this  is,  in  effect,  implied  by  the  WKB  representation  which  requires 
that  Sx/Sz  be  very  small.  In  fact,  if  particles  are  to  be  reflected 
frcm  the  surface  of  the  stream  (Ferraro  showed  that  this  is  to  be  ex¬ 
pected),  then  the  density  must  become  infinite  at  the  surface  in  order 
to  satisfy  the  continviity  relations.  This  will  be  discussed  in  greater 
detail  in  IV. A. 

The  net  effect  of  this  Implied  assumption  of  constant  density  is 
that  Ferreuro  obtains  the  boundary  condition  at  z  =  z^  : 


H^/8«  =  n  m^W^  ,  (28) 

which  is  not  a  true  momentum  balance.  Ferraro  e]q>lains  the  discrepancy 
thus  (p.  36):  '^Ince  the  motion  of  the  particles  Is  reversed. . .the  change 
of  momentum  is  2x01^..,.  The  result  Implies  that  one-half  of  the  momen¬ 
tum  of  the  stream  is  taken  up  by  the  magnetic  field,  and  the  other  half 
by  the  electrostatic  field. "  Since  the  electrostatic  field  is  internal 
to  the  pleisma,  this  is  not  a  satisfactory  e3q>lanatlon.  Dungey  has  shown 
that  the  proper  relation  is  Indeed 
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} 


(29) 


H^/8jt  =  2rm^ 


and  the  difference  is  due  to  the  implied  assumption  of  constant  density. 
Ferraro  also  inferred  that  for  an  oblique  flow 

H^/8«  =  n  cos^  X  ,  (30) 

when  the  flow  is  inclined  at  an  angle  X  to  the  z-axis.  (Note;  cos  X 
is  not  squared  in  Ferraro's  Eq.  (159) J  this  is  a  typographical  error, 
and  it  is  corrected  in  Ferraro,  1960b). 

It  should  be  noted  that,  for  constsuit  X, 

A  =  Aq  exp  [-  X(z  -  Zq)]  ,  (31) 

a  typical  electromagnetic  shielding  formula,  so  that  l/X  gives  a  lower 
(no  dissipation)  estimate  for  the  "skin  depth"  or  thickness  of  the  inter¬ 
face  between  the  GMF  and  the  SCR.  This  is  in  agreement  with  the  later 
work  of  Dungey  on  the  steady-state  case, 

B.  MAIN  PRASE  OF  STORM 


1.  Ring  Currents 

The  work  described  so  far  has  served  only  to'  establish  theoretical 
foundations  to  explain  the  augmentation  of  the  GMF  during  the  first 
phase  of  a  geomagnetic  storm.  However,  this  does  not  explain  the 
observed  decrease  in  the  GMF  (main  jdiase)  which  follows  within  a  few 
hours  of  sudden  commencement  (SC).  To  account  for  this  Chapman  and 
FerrarO;  (1933),-  postulated  a  westward- flowing  equatorial  ring  current 
at  several  earth  radii  due  to  charged  particles  penetrating  the  GMF 
from  the  approaching  SCB.  The  effect  is  to  produce  a  magnetic  field 
opposed  to  the  GMF;  however,  among  other  things,  the  theory  does  not 
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explain  the  time  lag  or  the  requisite  stability  of  the  circular  ring  in 
an  assumed  magnetic  dipole  field.  Nevertheless,  it  did  provide  a  reason¬ 
able  model  for  the  main  phase  decrease. 

In  195T>  Singer,  in  a  brief,  but  comprehensive^ review  of  magnetic 
storm  theories,  described  the  theory  of  the  interplanetary  shock  wave  due 
to  gas  erupting  from  the  sun.  The  impact  of  this  shock  wave  upon  the 
GMF  compresses  the  field  and  also  causes  a  magnetohydrodynamic  shock 
within  the  magnetosphere  which  registers  at  the  earth's  surface  at  the 
auroral  zones,  being  "channeled"  there  by  the  magnetic  field  lines.  As 
the  lines  of  force  converge  toward  the  poles  the  shock  wave  converges 
with  them,  increasing  in  strength  about  five  times  to  a  point  where  it 
can  produce  the  auroral  effects  generally  associated  with  SC.  Several 
hours  later,  the  gas  from  the  solar  flare  which  acts  as  the  "driver"  for 
the  shock  arrives  at  the  GMF,  penetrates  into  the  field  and  is  trapped 
in  the  magnetic  tubes  of  force.  What  follows  is  a  simple  application  of 
the  theory  of  charged-particle  drifts  in  anisotropic  magnetic  fields,  where 
the  anisotropy  may  be  Introduced  by  means  of  eui  electric  or  gravitationeQ. 
force,  or  through  the  inhcmogeneity  and  curvature  of  the  magnetic  field 
itself. 

Singer,  (1957);^  calculates  the  drift  of  protons  of  velocity  2000  km/sec 
injected  perpendicular  to  the  field  lines,  based  on  the  inhomogeneity  of 
the  dipole  field,  neglecting  magnetic  effects  due  to  the  peurticle  drifts 
themselves.  Sy  comhlnlng  this  motion  with  oscillation  along  field  lines 
he  finds  a  westward-flowing  ring  current  of  roughly  diamond-shaped  cross- 
section,  with  maximum  density  in  the  equatorial  plane  and  6u:ound  8R, 
extending  30°  north  and  south  of  the  equator.  Despite  the  rou£^ess  of 
the  csLlculations  they  clearly  show  that  particle  drifts  can  give  rise  to 
magnetic  fields  sufficient  to  produce  the  main  phase  of  the  storm. 

There  still  remains  doubt  as  to  the  cause  of  the  drifts:  how  can  the 
SCR  penetrate  the  CA(F  if  the  latter  has  already  been  coiiq>re8sed  due  to 
the  prior  arrival  of  the  shock  wave?  Singer  only  says  that  particles 
may  leak  into  the  field  if  it  is  sufficiently  distorted  from  the  magnetic 
dipole,  presvimably  in  hi^er  latitudes  where  the  CMF  disturbance  depends 
more  strongly  on  loced  time.  Dessler,  et.  eO..  I96I;  removed  this  objection 
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maintaining  that  penetration  is  not  necessary:  instead,  large-amplitude 
hydromagnetic  waves  are  generated  by  the  impact  of  the  SCR  upon  the  GMP; 
the  resulting  impulse  is  transmitted  in  the  form  of  longitudinal  hydro- 
magnetic  shock  waves  traveling  toward  the  earth,  and  treuisverse  waves 
traveling  toward  the  poles  silong  magnetic  field  lines. 

Like  ordinary  waves  in  a  compressible  fluid  (Liepmann  euid  Roshko, 

1957,  Chap.  3)^  these  shock  waves  tend  to  steepen  as  they  propeigate,  un¬ 
til  the  extreme  temperature  and  density  gradients  across  the  shock  dis¬ 
sipate  energy  fast  enovigh  to  prevent  further  steepening.  Thus,  as  the 
shock  passes  through  the  gas,  the  gas  is  heated,  and  the  higher  thermal 
energies  of  the  gas  particles  then  give  rise  to  higher  drift  velocities 
in  the  GMF,  due  to  field  inhomogeneities  and  cvirvature  of  field  lines. 

The  latter  also  gives  rise  to  a  westward-flowing  ring  current.  The 
heating  goes  on  until  the  gas  pressure  is  comparable  to  that  of  the  mag¬ 
netic  field  perturbation,  which  in  tvim  is  eciuEJ.  to  the  change  of  momen¬ 
tum  of  the  incident  SCR.  Thus  the  gas  particles  actjuire  thermal  velocities 
of  the  order  of  the  velocity  of  the  SCR;  the  resultant  situation  is  similar 
to  that  which  Singer  obtains  by  means  of  injection. 

When  the  ratio,  6B/b,  of  pertxirbed  field  to  steady  field  becomes 
much  less  than  vinlty,  the  hydromagnetlc  heating  decreases  rapidly,  so 
that  the  authors  estimate  the  maximum  ring  current  density  to  occ\u:  around 
4r  .  The  c8Llculatlons  which  substantiate  the  theory  are  only  ordei<-of- 
magnitude,  and  the  theory  offers  ample  scope  for  further  research. 

Although  the  existence  of  a  ring  current  to  the  source  of  the  main 
phase  decrease  is  still  in  doubt,  Slutz,  (1962).  has  noted  that  if  it 
does  exist  the  pressure  of  its  magnetic  field  would  enlarge  the  size  of  the 
geomagnetic  cavity,  since  it  would  augment  the  GMF  outside  the  ring.  Thus, 
we  find  the  minimum  dimension  of  the  cavity  (subsolar  point,  equatorial 
plane)  from  the  boundary  condition  [^.  (29)],  where  H  equals  the  dipole 
field  plus  the  Induced  field  due  to  surface  cwrents,  which  we  sheLLl  take 
equal  to  the  original  field 

(R/r^)^  =  (8«)  2nm^W^  ,  (31) 
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where  :  =  0.3  gauss.  This  yields 

Tq  =  8.85R  ,  (32) 

■which  is  smaller  than  the  experimental  estimates  which  are  in  the  neighbor¬ 
hood  of  13. 6r  .  However,  if,  on  the  basis  of  the  rather  fragmentary 
Explorer  VI  and  Pioneer  V  results,  we  assume  a  toroidal  ring  current  in 
the  equatorial  plane  and  coaxial  with  the  GMF,  with  radius  9R  SJid  cross- 
section  radius,  3R  >  and  cxirrent  of  5  x  10^  amps.  ,  then  Eq.  (31)  is 
altered: 

4[HQ(R/rQ)3  +  =  (8«)  2nmy  ,  (33) 

again  assuming  the  effect  of  the  sixrface  currents  is  to  double  the  original 
field  at  the  boxindary.  If  we  calculate  from  the  formula  for  a  ciar- 

rent  loop  we  find  that  r^  =  12.2  R  (S^^ythe,  1950,  p.  27I). 

2.  Aurorae 

We  have  seen  in  the  preceding  section  how  Singer  and  Dessler  could 
explain  auroral  activity  in  terms  of  the  propagation  of  hydromagnetic 
shocks  8G.ong  field  lines  into  the  auroral  regions.  One  need  not  have  a 
gecxnagnetic  storm  to  accoc^lish  this;  it  would  be  sufficient;  for  the  inter¬ 
face  to  be  easily  distorted  under  the  influence  of  random  fluctuations 
in  the  steady-state  SCR.  Thus,  fluctuations  above  some  minimum  thres¬ 
hold  energy  might  generate  the  hydromagnetic  waves  which  give  rise  to  the 
"steady-state"  aurora  which  occurs  most  of  the  time. 

It  has  also  been  argued  that  'there  are  "neutred  points"  in  speuze 
corresponding  to  -the  points  "N"  in  Fig.  3^  at  which  SCR  particles  can 
enter  the  magnetosphere  and  spiral  down  the  field  lines  toward  the  poles, 
where  they  give  rise  to  the  aurora  by  collisions  with  atmospheric  con¬ 
stituents.  This  does  not  explain  the  fact  that  auroral  activity  is  a 
maximum  on  the  nlghtslde,  \inle88  the  fields  are  sufficiently  twisted 
neeu*  the  poles,  due  to  the  rotation  of  the  earth. 
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Chapman;,  (1961a);  believes  that  high-energy  particles  in  the  SCR, 
especially  during  storms,  may  break  throvigh  the  GMF-SCR  interface  and 
become  trapped,  replenishing  and  distorting  the  Van  Allen  belts.  Then, 
at  those  points  on  the  night  side  of  the  magnetosphere  where  the  ring 
current  field  ceuicels  the  (MF,  one  finds  internal  neutral  points  where 
psurbicles  from  the  radiation  belts  may  spiral  down  to  the  earth.  These 
neutral  points  form  long  luDesbaped.  strips,  which  explains  why  the  a\irora 
shows  a  characteristic  ribbon-like  structiure.  The  dimensions  of  these 
strips  may  be  100  x  1000  x  1  miles  .  On  the  day  side  there  are  no  neu¬ 
tral  points  or  strips  since  the  GMP  is  con^ressed  by  the  SCR. 

There  is  as  yet  no  generally  accepted  theory;  the  question  remains 

open. 
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IV.  STEADY-STATE  GMF-SCR  INTERACTION 


A.  ONE-DIMENSIONAL  MODEL 

(1958^; Chap.  8).  applied  Ferraro's  plane  model  to  the  cal¬ 
culation  of  the  steady-state  interaction,  in  which  h/ht  =  0  and  the 

plasma-field  interface  is  stationary  in  time,  i.e.,  dz^/dt  =  0  .  We 

shall  here  expand  Dungey's  analysis  in  detail,  since  it  constitutes  bur 
best  mathematical  Justification  for  some  basic  ideas  about  the  inter¬ 
action  between  the  GMF  and  the  SCR.  In  addition,  we  present  calcialations 
of  particle  trajectories  and  currents. 

Following  Ferraro's  model,  consider  a  semi-infinite,  neutral,  fully- 
ionized  stream  of  protons  and  electrons  moving  in  the  negative  z-direction 
as  shown  (Fig.  5).  The  particles  are  taken  to  be  at  absolute  zero  temper¬ 
ature,  and  penetrate  into  the  permanent  magnetic  field  as  far  as  the 
plane  surface  of  the  plasma  at  z  =  z^  . 

We  asshme  that  points  in  the  y-direction,  and  is  given  by 

Hp  =  H^Ca/z)^  ,  (34) 

where  Hq  =  0.3  gauss  which  is  augmented  by  the  magnetic  field  H' 
in  the  y-direction,  due  to  surface  currents  in  the  z  =  plane. 

We  assume  the  existence  of  stationary  electric  and  magnetic  fields 
in  the  plasma,  which  may  be  represented  by  the  potentials 

♦  -  ♦(z)  ;  A  -  A(z)  ,  (35) 

where  A  is  the  vector  potentleuL,  and  is  chosen  to  be  pointing  in  the  x- 
dlrection,  so  that  the  magnetic  field  H  *  curl  A  within  the  plasma 
will  be  parallel  to  +  H'  ,  and  boundary  conditions  can  be  satisfied. 

The  particles  will  be  moving  at  z  ■  »  with  only  a  z-conrponent  of 
velocity;  there  eure  no  forces  (E  =  b^/bz  is  in  z-dlrection)  which  can 
produce  a  y-con^Kjnent,  so  that  the  velocity  of  a  particle  is  represented  by 
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V  * 


(u,  0,  w)  . 


(36) 


If  q.  Is  the  charge  on  a  particle  In  the  plasma  and  m  Its  mass,  Its  Hamiltonian 
is 


^  =  (2m) (p  -  +  (!♦  =  ^  (u^  +  w^)  +  q«t»  ,  (37) 

and  since  does  not  depend  explicitly  on  time  its  veilue  is  a  constant; 
thus 


'H  =  (1^2)  (U^  +  W^)  ,  (38) 

vhere  W  is  the  particle  velocity  in  the  negative  z-direction  and  U  is  the 
velocity  in  the  x-direction  at  z  =  »  ,  the  same  for  electrons  and 

protons. 

Since  the  densities  of  Interest  are  so  low,  one  is  Justified  in 

neglecting  inter-particle  collisions;  hy  also  neglecting  any  radiation 

throvigh  hremstrahlimg  we  have  constructed  a  conservative  Hamiltonian  in 
2  2 

which  (u  +  w  )  depends  only  on  z.  The  property  of  tlme-reVerslblllty 
shows  that  the  particle  trsijectories  will  be  symmetric  with  respect  to 
a  line  peorallel  to  the  z-axis  drawn  throu^^  their  tuming;<poists  in  the 
Zq  plane.  Thus  the  absolute  veilue  of  w(z)  for  both  incoming  and  out¬ 
going  particles  of  a  single  type  is  the  same,  althou£^  different  in  sign; 
however,  u(z)  does  not  change  sign. 

The  charge  density  at  some  value  of  z  is  due  to  the  density  of  both 
incoming  and  outgoing  particles.  For  charge  density  p  and  current  density  J: 

P  =  2e(n  -  n  )  ;  J  »  |S;(n  u  -  n  u  )  ;  J 

le  X  c  .iiee  y 

Maxwell's  equations ' take  the  form 

dlv  E  *  4«p  >»  8«e(nj^  -  n^)  ,  (4o) 


(39) 
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and 


curl  H  =  4jtJ  =  8jt  —  (n.u.  -  n  u  )  .  (4l)  / 

X  c  '  i  i  e  e'  / 

t 

In  the  gecmagnetic  region,  z  <  ,  the  existence  of  a  permanent 

magnetic  field  can  be  thought  of  as  due  to  a  distribution  of  permanent 
magnetization,  or  cxirrent  dipoles,  which  has  the  vector  potential  (in  the 
x-direction) : 


Ap  «  "  i  Hq  a(a/z)^  .  (42) 

The  field  is  aiagmented  by  H'due  to  the  current  in  the  plasma.  Since 
there  are  not  so\irces  of  Il'in  the  geomagnetic  region,  and  since  symmetry 
requires  that  its  vector  potential  be  in  the  direction  of  x  euad  dependent 
upon  z,  the  magnetostatic  equation 

^  =  0  ,  (43) 

dz 


yields 

A'  «  H'z  , 
H'  *  const. 


(44) 


Finally,  when  the  solution  is  derived  the  following  boundary  conditions 
must  be  satisfied: 

(i)  vector  potentiaG.  continuous  across  the  boundary 
z  -  Zq  ; 

(li)  magnetic  field  intensity  continuous  across  the 
boiindeuy; 
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(iii) 

(iv) 

(v) 

(vl) 


electrostatic  potential  continuous  across  the 
boundary j 

conservation  of  energy; 

over- all  charge  neutrality; 

no  fields  at  z  =  »,  where  ion  and  electron 

densities  are  equal  to  n^,  and  velocity  equals 

W  in  the  negative  z-direction,  U  in  the  x-directlon  . 


The  unknown  quantities  are  the  velocities,  potentials,  and  densities; 
they  are  first  derived  ae  functions  of  the  potential.  A,  beginning  with 
the  lateral  velocities.  Hamilton's  equations  of  motion  yield 


0 


> 


whence 


p^  =  mu  +  qA  =  const, 
and  condition  (vi)  requires  that 


(J^5) 


U  -  u  =  +  ^  ,  (46) 

me  '  '  ' 

which  is  position-dependent,  as  assumed  in  Eq.  (35)' 

The  velocity  cooq>onents  in  the  z-direction  can  be  found  from  the 
original  Hamiltonian,  Eq.  (37): 

(mj^/2)w^  =  (mj^/2)W^  +  eAU/c  -  e^A^/2m^c^  -  e#  (4T) 

(m^2)w^  *  (m g/2)W^  -  eAO/c  -  e^A^/2m^c^  +  e#  ,  (48) 

where  we  have  substituted  for  u  from  Eq.  (46).  We  now  make  the  assun^stion 
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of  "quasi -neutral! tyj"  i.e.,  an  electrostatic  field  due  to  charge  separa¬ 
tion  does  exist,  represented  by  -  grad  ♦,  but  a  sufficiently  strong 
field  may  be  obtained  from  relatively  little  change  in  density,  i.e., 

n^^  =  n^  =  n  .  (49) 

In  the  steady  state  condition  the  divergence  of  the  electric  current 
density  must  be  zero  in  both  inccming  and  outgoing  streams,  so  that 
Eq,.  (49)  and  condition  (vi)  lead  to 

=  w  .  (50) 

Since  Eq,.  (49)  would  appear  to  violate  Eq.  (4o),  the  criterion  for  validity 
of  this  approximation  is 


div  E  . .  T 


(51) 


The  results  of  the  calculation 
If  we  substitute  Eq.  (50) 


are  consistent  with  this  criterion, 
into  Eqs.  (47)  and  (48),  we  find; 


♦ 


A  .i^U^U/c 

..2  m_  m^  ' 


(52) 


and 


(53) 


2  .2  2/2  2 
where  ■  m^^m^w  c  ye  Is  Just  the  value  of  A  when  z  ■  z^  , 

derived  by  setting  w  ■  0  in  Eqs.  (47)or  (48).  For  simplicity,  we 

shall  set 
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a  s 


a(z)  }  Aq  <  0  ; 


ot  will  prove  to  be  the  most  convenient  V8u:lable  for  this  analysis. 

The  particle  density  is  foiind  from  the  equation  of  continuity 
which  must  hold  for  each  stream  separately,  so  that 

nw  =  n^W  ;  (5^) 


therefore, 


^  .  (1  -  .  .  (55) 

0 

It  is  now  possible  to  calcvilate  o(z)  by  substituting  Eqs.  (^*6)  and  (55) 
into  Eq.  (4l).  The  result  is 


curl  H  * 


(56) 


where  M  is  the  reduced  mass,  M  *  m^m nig)  •  The  first  integra¬ 
tion  of  this  equation  subject  to  condition  (vi)  yields 


tl  -  (1  -  a^)^)  • 


(57) 


Since  H  »  ,  the  condition  which  determines  the  boundary 
(ot  *  l)  is  that  incoming  particles  are  specularly  reflected  by  the 
pressure  of  the  total  magnetic  field  (see  Fig.  5): 


(58) 


2 

“e^  ^  ^  ^ * 

It  should  be  noted  that  Ferraro  finds  H^/8rt  =  1  ^  ^  * 

which  violates  conservation  of  momentum.  This  is  beca\ise  Ferraro  in 
essence  considers  n  =  constant,  and  does  not  account  for  its  sudden 
increase  at  the  boundary;  this  assumption  leads  to 


X.^  =  SjoiQe^/Mc^  ;  d^Ot/dz^ 


;  a  =  e’^^^  '  ,  (59) 


where  l/X  is  the  ion  Larmor  radius  at  z  =  . 

The  fields  in  the  region  z  <  z^  may  now  be  determined  by  conditions 
(i)  and  (ii),  which  yield 

•  &l  Itl’  ■  ^ 

H'  +  j-l-l  =  4[w  nQ(mj^  +  .  (6l) 


Since  z^  >  a  ,  a  sin^jle  numerical  cedculation  shows  that  the  right 
side  of  Eq.  (6o)  is  negligible  ccsnpared  with  the  magnitvides  of  the  terms 
on  the  left,  resulting  in 


1^1  ^  ■ 
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To  find  the  vector  potential  in  the  region  z  >  z^,  Eq.  (57)  is  integrated: 

(e^)  exp  [-X  (z  -  Zq)]  =  tern  '^o®  (0/2)]  ,  (64) 

where  sin  6  =  cc  .  Figure  6  compares  Perretro's  exponential  solution 

with  Dungey's  exact  solution. 

Explorer  X  ( IGY  Bull . ,  1962a,  p.  13 )  observed  a  flux  of  SCR  of 

n  =  15/cm^  and  W  =  7OO  kn/sec  which  coincided  with  the  sudden 

commencement  of  an^lltude  3OO7  observed  at  College,  Aleiska.  Under 
these  conditions  Eqs.  (62)  and  (63)  yield  H'  =  8O7  when  the  front 
of  the  storm  stream  comes  to  rest  in  the  GMF.  This  constitutes  a  fair 
approximation,  considering  that  the  actiuG.  field  is  curved,  not  planar, 
and  would  be  magnified  in  the  aviroral  zones  according  to  Singer  (see 
III.  B.l). 

If  we  consider  the  trajectory  of  a  single  incoming  particle  of 
charge  q,  mass  m,  with  time  coordinate  -oo<t<»  in  a  system  moving 
with  V  »  (U,  0,  0)  ,  then  the  arc  in  space  is  described  by 


Equation  (56)  allows  us  to  replace  the  right. side  of  Eq.  (65)  by 

dx  /  IAq  \ 
dz  ImxSfc  dz^ 

^mX'Tfc/  \dzj 

and  substituting  into  (67)  from  Eq.  (57).-  yields 


(66) 

(67) 
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PIG.  6 — Comparison  of  exact  solution  for  Of  with  exponentleil  approximation 
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} 


(68) 


X,  = 


m  X 
e  e 


(69) 


The  outgoing  trajectory  is  syinmetrical  with  respect  to  the  turning  point 
(see  Pig.  7). 

From  Eq.  (56),  for  the  current  density  ,  ve  have 


a(i  -  c?)"^  >0 


(TO) 


and  the  total  current  per  xmit  length  of  plasma  front,  measured  in  the 
y-direction  is 


2e 

He 


j 


«(1 


)^  dz 


(71) 


which  is  most  easily  integrated  in  terms  of  0  <  ot  <  1  ,  hy  substituting 

for  dz/dP!  from  Eq.  (56).  This  gives 


[l6«nQ 


(72) 


nine-tenths  of  which  flows  near  the  surface,  in  the  region  where 
0.l4l  <  ®  <  1.  .  This  is  identical  to  the  eaqpression  for  surface  cur¬ 
rents,  if  one  were  to  assume  con^lete  sepeuration  of  magnetic  field  and 
plasma. 
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FIG,  7- -Electron  trajectory 
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In  order  to  verify  the  approximation  of  quasi -neutrality,  we  calculate 
the  charge  density  from  Poisson's  Equation,  which  hecomes 

-  d^V^z^  =  4rtp  i  -  (ejm^eW^/Mc^)  [(m^  -  m^)  (2  -  ^  1  -  +  a^/Vl  - 

+  (m^m^)^  (U/W)  a/Vl  -  .  (73) 

When  this  is  integrated  over  z^  <  z  <  »  hy  means  of  the  substitution 
dz  =  (dz/dO!)  dot  ^  we  find  a  net  charge  in  the  plasma  per  imit  surface 
area  of 


Q  =  -  (m^  c^Uq/  nm^)^  (w/c)^  [l  +  (u/W)  (m^m^^)^]  ,  (7^) 

approximately,  to  within  (m ^m^  «  l). 

To  preserve  overall  neutredity,  we  require  a  positive  surface  charge 
Of  (-  Q)  at  z  =  Zq  .  This  meeuis  that  all  electric  lines  of  force 
which  originate  in  the  plasma  end  on  the  positive  charges  in  the  surface, 
so  that  there  is  no  electric  field  in  the  region  z  <  ,  and  the  elec¬ 

tric  potential  there  is  a  constant.  The  surface  charge  coiresponds  to 
physical  reality  in  that  the  more  massive  ions  will  penetrate  further 
into  the  magnetic  field  than  the  electrons,  while  being  retarded  by  the 
Coulomb  force  because  of  their  sll£^t  sepeuratlon  from  the  electrons.  We 
may  verify  the  quasi -neutrality; 


neglecting  the  term  in  (U/W)  gives 


<  0.005 


(76) 
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at  the  surface  of  the  plasma,  and  generally  much  smaller  than  O.OO5  in 
its  interior.  ■ 

A  more  general  discussion  of  this  problem,  Including  consideration  of 
(one -dimensional)  velocity  distributions  at  »  =  x  is  to  be  found  in 
Grad's  I96I  paper. 

B.  TWO-DIMENSIONAL  MODEL 


We  now  formvilate  an  idealized  problem,  using  the  insights  gained  in 

section  IV. A.  We  consider  the  two-dimensionsQ.  problem  first,  i.e.,  a 

configuration  of  line  currents,  located  at  the  origin  of  coordinates,  and 

small  in  extent  compared  to  the  dimensions  of  the  problem.  We  wish  to 

find  the  magnetic  field  within  a  certain  region  containing  the  origin;; 

the  bovindary  of  this  region  is  to  represent  the  interface  between  the  GMF 

and  the  SCR.  Within  the  region  both  the  scalar  (G)  and  vector  (A  )  mag- 

z 

netic  potentials  must  satisfy  Laplace's  equation: 

with  a  specified  singrilarity  at  the  origin  due  to  the  given  configuration 
of  line  current  sources. 

In  the  preceding  section  we  found  that  the  magnetic  field  penetrates 
only  a  very  short  distance  into  the  plasma  stream;  in  our  idealized  model 
we  assume  that  the  field  is  completely  excluded  from  the  plasma  by  c\ir- 
rents,  J,per  unit  length,  flowing  in  the  Interface.  These  currents  must 
satisfy  the  relation 

H  -  4«J  ,  (78) 

where  H  is  the  total  field  at  the  bo\andary.  Since  normal  coaq^onents  are 
continuous  across  the  boundary  we  have  the  first  boundaiy  condition: 

FpA-C*0;  C*  const.  , 
z 
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vhich  is  equivalent  to  the  requirement  that  the  magnetic  field  he  parallel 
to  the  interface.  If  C  can  he  determined  at  one  point  on  the  boundary 
hy  considerations  of  symmetry,  for  example,  and  if  the  hovindary  vere  knovn, 
ve  should  have  a  Dirichlet  problem,  which  has  a  unique  solution. 

The  boundary  is  not  known;  instead  we  have  the  second  condition 

cos^  X  ;  =  8jt(2nQV^)  +  nig)  >  P  >  0  ;  (79) 

where  X  is  the  angle  between  the  normEti  to  the  boundary  and  the  x-cucis 
which  is  now  taken  to  be  the  direction  pareOlel  to  the  stream. 

This  second  boundary  condition  must  compensate  for  our  lack  of  know¬ 
ledge  concerning  the  exact  location  of  the  bounding  surface.  We  shall 
assume  that  the  problem  thus  formulated  is  susceptible  of  a  unique  solu¬ 
tion;  however,  this  has  not  yet  been  rigorously  established,  to  my  know¬ 
ledge,  and  would  constitute  a  valuable  Eiddition  to  the  theory  of  free¬ 
boundary  problems. 

2  2  2  2  2 
Ve  may  express  cos  X  =  (dy/ds)  ,  where  ds  »  dx  +  dy 

is  the  arc  length  along  the  boundary  in  the  xy-plane  (see  Fig.  9).  Thus, 

the  boundary  condition  mey  be  expressed  as 

-dfl  »  Hds  =  ±Pdy  ,  (80) 

with  i  depending  upon  whether  the  field  lines  sure  parallel  or  antipeurallel 
to  ds  .  Over  some  particular  region  of  the  Interface  we  may  then  Integrate 
the  above  equation: 


0  .  a±py-K  *  0  (8l) 

K  B  const. 

Thus,  the  Integral  form  of  the  two  boundary  conditions  defines  two  one- 
parameter  families  of  surfaces,  F  and  0  ,  the  intersection  of  which  is 
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-The  tvo-diaenslonal  line  cvirrent 


the  desired  interface.  In  two  dimensions  and  A  are  Cauchy -Riemann 

z 

conjugates,  so  that  the  choice  of  a  representation  for  one  automatically 

determines  the  other.  One  method  of  obtaining  an  approximate  solution 

(see  below)  would  be  to  represent  0  and  A  by  an  infinite  harmonic  series 

z 

with  \indetermined  coefficients  pl;is  separate  terms  to  represent  the  source 
of  the  field.  One  could  then  evaluate  the  coefficients  so  that  the  two 
surfaces  would  be  made  to  fit  as  closely  sis  desired  over  the  region  of 
interest  (co-location  method).  Again,  it  must  be  noted  that  the  validity 
of  the  foregoing  has  not  been  rigorously  determined. 

Both  Dungey  (196I)  and  Zhigulev  ( 1959c)  have  attempted  to  solve 
the  problem  of  a  two-dimensionsd.  magnetic  line  dipole  in  the  streaming 
SCR,  with  the  aid  of  various  conformal  transformations.  However,  the 
most  lucid  derivation  is  that  of  Hurley  (l96la,b)  who  has  treated  both 
the  problem  of  a  line  c\irrent  transverse  to  the  direction  of  the  SCR, 
and  the  problem  of  the  line  dipole  inclined  at  surbitrsuy  angles  to  the  SCR. 

In  the  case  of  the  line  cxirrent,  I,  Hurley  shows  that  the  equation 
of  the  interface  is 


2  cos  (Py/^l)  =  exp(Px/4l)  ;  (82) 

the  magnetic  field  may  be  derived  from  either  the  magnetostatic  scslIsu: 
potential 

n  =  l(y'  +  surcsin  U)  ,  x'  =  Px/Ul  (83) 

U  =  (1  -  cosh  x'cos  y')/(co8  y'  -  cosh  x')  ,  (84) 

or  the  vector  potential 

A  =  -  l[x'  +  log(2  cosh  x'  -  2  cos  y')]  .  (85) 

z 

In  Fig.  9  we  have  normalized  to  the  smallest  dimension, 

^0  “  loSgS  .  (86) 
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.  By  meems  of  the  co-location  method  descrihed  above  one  can  obtain 
the  series  representation  of  the  F  and  G  surfaces  to  second-order  approxi¬ 
mations: 


F  =  log^Cr/r^)  +  0.3^^73(r/rQ)  cos  6  +  0.01993(r/i’Q)^  cos  26  -  0.36466  = 

(87) 

G  i  6  -  1.03986(r/rQ)  sin  6  +  0.01993(r/rQ)^  sin  26  =  0  (88) 

where 

r^  =  2.7721/P  .  (89) 

Fig\ire  9  shows  the  F  and  G  surfaces  as  well  as  Hurley's  exact  sol¬ 
ution.  Magnetic  field  lines  inside  the  Interface  have  also  been  plotted. 
It  appears  that  substantial  agreement  between  the  F  and  G  surfaces  in¬ 
dicates  a  good  approximation  to  the  rijSht  answer;  it  may  also  prove  pos¬ 
sible  to  generate  a  quantitative  estimate  of  the  actual  error  based  upon 
the  difference  between  F  and  G. 

Hurley  finds  H'  =  O.^H^  ,  H  «=  H'  +  H^  ,  at  (r^^O)  which 

is  comparable  to  the  vedue  H'  *  derived  from  the  Ferraro-Dungey 

plane  model. 

The  case  of  a  line  dipole  inclined  at  angle  ct  to  the  y-axls  is  more 
difficult,  and  the  explicit  expression  for  the  fields  or  their  potentleG. 
is  not  derived  by  Hurley,  althou^  it  is  fully  Implicit  in  his  solution. 
Thus,  if 


H^  -  iH^  *  H*  e  d«/dz  ;  *  »  -  G  +  lA  ,  (90) 

X  ^  z 

then  the  nature  of  the  singularity  at  the  origin  is  such  that 

♦  ->2lViz  as  2  -►O  .  (91) 
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Hiirley  transforms  to  the  w-plane,  such  that  the  houndary  is  mapped  onto 
a  circle  of  radius  a;  i.e.,  v  =  a  e  •  ?y  requiring  that  •w(z)  -*-z 
as  z  -+0  he  derives 

4  =  (2M/i)  (l/w  -  w/a^)  ,  (92) 

for  the  field  of  a  dipole  boaonded  by  a  circle;  he  then  finds  the  trans¬ 
formation  z  =  z(w)  vhich  distorts  the  boundary  of  the  circle  until 
the  pressure  condition  is  satisfied.  Hurley  then  shows  that  for 

2 

a  =  ;  M  =  dipole  moment  ,  (93) 

and 

z'  =  zO/itM)^  , 

the  equation  of  the  boundary  is  given  in  terms  of  ♦  eus  parameter: 

lac'  =  -  [l  -  8in(4  +  a)]  ,  logtcos  (♦  +  o:)/(i  +  sin(4  +  ot))]^  (9^) 

+  log[(l  +  cos  ♦)/(!  +  8in(4  +  Ct))]^  +  2 

8in(4  +  a)  +  2  -  aa/n  for  «  >  ♦  >  rt/2  -  a  (95) 

y'  =  8ln(4  +  a)  -  act/*  for  jt/a  -  o  >  ♦  >  -  n/2  -  a 

-  8in(4  +  a)  -  2  -  aa/n  for  -  jt/a  -  a  >  ♦  >  -  n  . 

We  have  applied  the  co-location  method  to  the  case  of  a  *  n/2  , 

using  the  third-order  approximation  thus  obtedned  as  a  trial  surface  in 
a  relaxation  calcvilation.  To  apply  the  relaxation  method  one  considers 
the  Dlrlchlet  problem  for  the  vector  potential  with  a  given  trial  sur¬ 
face.  If  the  approximation  is  a  reasonably  good  ope,  only  relatively  snail 
perturbations  of  the  boundary  surface  are  required  to  obtain  agreement 
with  the  second  boundary  condition.  In  this  case  the  trial  surfaces 
quickly  converged  to  Hurley* s  solution. 
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C.  THREE-DIMENSIONAL  MODEL 


To  the  present  time  (July,  19^2)  there  have  "been  only  a  half- 
dozen  published  papers  •which  give  same  quantitative  estimate  of  the 
three-dimensional  interface  euround  a  stationery  magnetic  dipole  per- 
pendlculeur  to  the  flow  of  the  SCR.  The  first  to  deed,  with  this  subject 
was  Beard  (i960)  who  represented  the  interface  by  a  function  F(r,  0,  t) 
derived  by  eui  approximate  method  from  the  second  (pressure)  boundary 
condition  expressed  in  the  form  of  a  nonlinear  partial  differential 
equation  in  F.  The  first  (-tangential)  boundary  condition  is  Ignored 
entirely  and  replaced  by  the  assumption  that  the  original  dipole  field 
is  only  changed  in  magnitude  by  a  constant,  k,  at  the  boundary,  without 
any  change  in  direction;  the  component  of  the  field  normal  to  the  boundary 
is  sin^ily  Ignored.  In  fact,  this  is  a  fruitful  source  of  error  since 
the  curvature  of  the  boundary  necessarily  results  in  considerable  dis¬ 
tortion,  aa  well  aa  augmentation  of  the  GMF.  However,  provided  one 
assumes  the  appropriate  value  of  k  (Beard  takes  k  >  2)  this  does 
constitute  a  surprisingly  good  approximation,  at  least  in  the  two-dimen- 
tional  case.  Beard  notes  that  near  the  poles  one  expects  to  find  neutral 
points  corresponding  to  Chapman's  N  (see  Fig.  3)  at  which  the  direction 
of  stream  flow  is  locally  parallel  to  the  surface,  so  that  H  b  0  ; 

he  eJLso  believes  that  in  this  region  the  ciuyature  of  the  surface  will 
be  much  higher  than  elsewhere,  and  may  even  require  modified  boundary 
conditions,  since  -these  are  based  on  the  eissuaptlon  that  the  ion  Larmor 
radius  is  small  relative  to  the  curvature  of  the  surface. 

Althou£pi  the  two-dimensloned  work  indicates  -the  usefulness  of 
Beard'*  s  first  supposition  it  Is  unlikely  -that  the  radius  of  curvature 
would  ever  approach  the  Lamor  radius;  if  It  did  one  would  have  almost 
a  cusp  In  -the  surface,  ideally  placed  to  trap  incoming  pauticles,  whose 
pressure  would  then  broaden  it.  Such  a  region  would  be  eig^ected  to  be 
highly  unstable  in  any  case  (Northrop  and  Teller,  I96O;  Chang,  I962). 
Beard  derives  an  expression  for  -the  shape  of  the  surface  in  the  ecurth- 
sun  plane  (dF/d4  s  0)  which  he  finds  to  be  approximately  circular. 

On  the  nlc^tslde  of  the  surface.  Beard  estimates  the  length  of  -the  tall 
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(see  Fig.  l)  to  be  approximately 


^  =  V/\  ’  (9^) 

where  is  the  thermal  velocity,  and  is  the  maximum  cross-section 
of  the  cavity,  the  tail  being  assumed  to  be  approximately  conical  in 
shape.  With  an  estimate  of  R-  <  1.8r^  ,  r^  the  minimum  dimension 

of  the  cavity,  approximately  r^  =  IQR  ,  and  T  =  10^  °K  ,  the 

tail  would  be  a  cone  of  eJ.titude  approximately  lOOR  . 

Sprelter  and  Briggs  (1962)  claim  to  have  discovered  inaccvuracies 
in  Beard’^s  work  and  have  re-done  it  obtaining  solutions  in  the  equatorial 
and  meridian  plane.  Unfortunately,  being  laiaware  of  Dungey'’s  work  at 
the  time,  the  authors  used  Ferraro's  bovindary  condition 

H^/8n  =  nmV^  . 

This  error  was  noted  and  rectified  in  a  later  letter,  indicating  that  their 
original  dimensionless  lengths  must  be  multiplied  by  (2)”^^^  when 
*  restored  to  physical  dimensions. 

More  recently  Beard  (1962)  has  published  euiother  paper  on  this 
subject  in  which  he  uses  the  first  approximation  to  calculate  surface 
currents  in  the  interface,  and  integrates  these  numerically  to  find 
their  effect  on  his  assvuned  dipole  field;  the  field  is  then  altered 
and  the  calculation  repeated,  thus  accounting  in  some  meetsure  for  the 
distortion  of  the  CMF  due  to  c\irvature  of  the  surface.  This  is  only 
calculated  for  the  meridlctn  plane,  which  once  again  yields  a  circular 
shape,  only  jf)  leurger  at  the  poles  than  at  the  equator. 

Slutz  (1962)  and  Midgeley  and  Davis  (1982)  have  treated  the 
case  of  a  dipole  enclosed  by  a  uniform  pressure  (  no  cos  X  in  second 
boxindary  condition).  Both  solutions  were  numerical,  but  by  different 
methods,  yet  the  results  for  the  cross-sections  agree  to  within  1^ 
at  the  pole,  and  3$  in  the  equatorial  plane. 
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D.  ANALOGIES  FROM  GASDYNAMICS 

A  conmon  technique  in  the  study  of  hypersonic  gaseous  flow  around 
solid  objects  Is  to  fix  the  object  In  a  wind  tunnel:  one  then  observes 
a  detached''st€uidlng”  shock  wave  on  the  upstream  side  of  the  test  object, 
the  geometry  and  thickness  of  the  shock  wave  depending  upon  Mach  number, 
temperatoire,  viscosity,  shape  of  object,  etc.  Downstream  of  the  object 
one  may  observe  a  "wake,"  an  area  of  relatively  low  pressure,  with  a 
sudden  change  In  fluid  velocity  at  its  boundary,  accompanied  by  turbulence. 

The  gas  In  passing  through  the  shock  front  Is  reduced  to  subsonic 
velocities,  regaining  speed  as  It  flows  around  the  object  and  Its  wake. 
There  Is  a  third  region  of  Importance,  and  that  la  the  boundary  layer 
between  the  object  and  the  gas.  .  This  is  a  transition  region 
between  a  region  where  no  gas  Is  allowed  to  flow,  l.e.,  the  object  Itself, 
and  the  region  of.  Irrotatlonal  flow;  it  Is  due  to  the  viscosity  of  the 
fluid  which  requires  that  It  be  stationary  at  the  surface  of  the  object. 
The  thickness  of  this  boundary  layer  increases  as  one  moves  downstream 
and  Is  proportional  to  the  square  root  of  viscosity  for  small  viscosities . 

The  GMF-SCR  interaction  has  been  compared  with  Just  such  a  problem, 
since  the  protons  are  hypersonic;  the  CWF  would  be  the  test  object,  and 
the  Interface  has  been  called  the  bovindary  layer.  On  the  basis  of  these 
analogies  It  has  been  suggested  that  there  should  be  another  transition 
region,  or  shock  wave.  In  the  SCR  between  the  eeurth  and  sun  due  to  the 
presence  of  the  geomagnetic  obstacle  to  the  flow,  and  also  a  "wake"  on 
ihe  night  side.  If  these  regions  do  exist  they  have  not  been  detected. 

There  Is  another  reason  for  doubting  the  existence  of  this  "detached" 
shock  wave:  If  we  calculate  values  of  density,  velocity,  ten^rature  and 
magnetic  field  upstream  of  the  presumed  shock  by  means  of  continuity  re¬ 
lations  using  values  measured  by  E;q;>lorer  X,  we  find  a  subsonic  to  super¬ 
sonic  transition  for  SCR  passing  through  the  shock.  Such  a  condition  Is 
forbidden  by  the  second  law  of  thermodynamics  In  the  case  of  both  clas¬ 
sical  (Llepmann  and  Roshko,  19^1,  pp.  36-61).  and  hydromagnetlc  shocks 
(de  Hofftnann  and  Teller,  1950#  PP«  698- TOO).  That  Is,  although  the  con¬ 
servation  of  energy  allows  subsonic  to  8iq>ersonlc  transitions,  the  re¬ 
quirement  of  Increasing  entropy  will  not  permit  Increasing  the  average 
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velocity  of  a  fluid  at  the  expense  of  its  thermal  energy  by  passing  it 
through  a  dissipative  region  like  a  shock  wave. 

Previously  we  have  considered  the  SCR  to  be  specularly  reflected  at 
the  interface;  every  measiarement  has  shown  considerable  turbxilence  at 
this  boundary,  and  it  may  be  more  correct  to  assume  diffuse  reflection 
at  the  boundary.  This  would  reduce  dependence  upon  the  angle  of  inci¬ 
dence,  but  not  greatly  alter  the  nature  of  the  problem.  A  more  extreme 
assun5)tlon  would  be  that  particles  can  actually  penetrate  the  transition 
region  which  separates  two  fairly  stable  regions,  the  magnetosphere  and 
the  SCR,  and  are  not  reflected  at  all.  Cahill  and  Amazeen  (1962)  find 
the  thickness  of  the  boundetry  to  be  100  to  1000  km  neeu:  the  subsolar  point, 
i.e.,  still  relatively  thin.  Under  this  assvmption  one  could  apply  the 
nonrelatlvlstlc  relations  of  de  Hofftaann  and  Teller  (1950)  for  a  hydro- 
magnetic  shock,  since  they  can  be  derived  by  Integrating  the  classic  con¬ 
tinuity  relations  across  the  transition  region  between  two  equilibrium 
states . 

For  a  normal,  one-dlmenslonELl  shock  in  a  transverse  magnetic  field, 
these  relations  are: 


mnu  =  A;Hu  =  B, 


(97) 


(a,  B,  C,  D  are  constants)  and 


p  +  mnu^  +  H^/8n  *  C  ;  E  +  u^/2  +  p/nm  +  1^/kiam  *  D  ,  (98) 


where  m  «  particle  mass,  E  >  internal  (thermal)  energy  per  milt  mass 
and  p  >  nkt  ,  pressure.  For  an  ideal  gas  we  also  know  that 
E  -  (3/2)nkT/  nm  ■  yi!Z/2aL  *  (3/2)p/nm  .  The  subscripts  (^)  or  (g) 

are  understood  since  the  above  equations  apply  on  either  side  of  the 
transition  region  (see  Fig.  lO). 
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Region 


PIG.  10 — Penetration  Model. 


The  constants  A,  B,  C,  D  can  he  evaluated  if  four  of  the  eight  un- 
knovn  quantities,  u,  n,  T,  H  are  specified.  In  this  case  let  us  assume 
the  quantities  upstream  of  the  transition  region,  or  interface,  to  he 
given  hy  the  Explorer  X  results: 


u^  *  300  lan/aec  ;  n^  »  lO/cm^  ;  ■  lOy  ; 

T^  w  0.5  X  10^  °K  . 

thus,  if  we  substitute  for  n,  T,  and  H,  the  last  equation  becomes  a  cubic 
in  u: 


4Au^-  5Cu^  +  2nAu  +  -  0  .  (99) 

By  symmetry  It  is  clear  that  u  «  u^  is  a  root  of  the  above  equatlcm; 
upon  dividing  by  (u  -  u^)  we  obtain  the  quadratic: 
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0 


(100) 


4  A  u  +  (4  Au^  -  5C)u  +  (2  DA  +  4 


Au^ 


50113^)  = 


of  which  only  the  positive  root  is  admissible, 

Ug  =  108  km/sec  , 


and,  therefore, 

=  27.8  j  =  27.87;  Tg  =  2  X  10^  ?K  ,  Pg  =  .76  x  lO"®  dy/cm^  , 

The  value  of  agrees  with  results  of  Explorer  X  and  Pioneer  I 
which  are  approximately  3O7.  However,  it  should  be  noted  that  the  data 
is  sparse  in  the  case  of  Pioneer  I  and  in  the  case  of  Explorer  X  we 
really  have  an  oblique  "shock"  at  some  \inknown  angle  to  the  direction 
of  the  flow.  Fxxrthermore,  in  preliminary  results  from  Explorer  XII  it 
was  mentioned  that  Hg  »  6O7  ,  presumably  measured  at  Euround  local 

noon.  The  value  of  was  not  given;  however,  it  was  stated  that  "the 
location  of  the  outer  edge  of  the  region  of  trapped  particles  appears  to 
be  at  or  below  the  termination  of  the  geomagnetic  field,  (IGY  Bull.,  1962b, 
p.  3)-  The  low-energy  proton  analyzer  with  a  threshold  of  200ey,  celling 
of  20  kev  ;  (IGX  Bull.,  1961b);  measured  no  protons  below  the  transition 
region.  However,  the  stream  energy  mu^2  ■  55  ev,  and  the:  thenoal-  energy 
corresponds  to  around  26O  ev  ,  so  that  measurements  neither  refute  nor 
substantiate  the  foregoing  ce0.culatlon.  Since  the  speed  of  a  compresslonal 
wave  propagating  normal  to  the  magnetic  field  is 

c^  =  5p/3nm  +  H^/4«nm  ,  (lOl) 


we  have 


a  1.035  X  10' 


2.02  X  10'  cn^sec 
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•with  Mach  nuaiberB,  respectively. 


Mj^  =  u^^/c^  *  2.9  ;  -  Ug/cg  »  0.535  • 

This  standing-shock  model  mi^t  explain  how  particles  penetrate  the 
magnetosphere,  euad  perhaps  account  for  ■the  experimental  measures  of  "the 
cavity  dimensions]  the  amount  of  mass  collected  by  a  cross-section  of 
ttCiOR)^  is  2  X  10  ear'bh  masses/lO^  years  .  If  the  GMF  penetrates 

all  space,  and  is  only  altered  near  the  interface  due  to  surface  currents, 
then  the  induced;,  field  H'  '  of  'these  currents  is  given  by 

*  (Hg  -  H^)/2  =  8.97  ,  (102) 

and 

GMF  =  (Hg  +  Hj^)/2  ;  GMF  »  l8.97  •  (103) 

The  earth’s  magnetic  dipole  field  will  have  -this  value  when  r  ■  12R  , 

which  is  closer  to  the  experimental  results  than  "the  Dungey-Ferraro 
estimate  of  8.1R  . 

It  should  be  noted  that  'this  concept  of  the  steady-state  standing 
shock  wave  is  different  from  •that  of  Singer  (1957)  who  believes  that  a 
bvu:st  of  SCR  lea'vlng  -the  sun  gives  rise  to  a  shock  wave  in  much  -the  same 
manner  that  a  "driver”  gas  causes  a  shock  wave  in  a  laboratory  shock  tube. 
It  is  this  traveling  shock  wave  which  causes  the  first  ihase  of  'the  geo¬ 
magnetic  storm  upon  striking  the  eeurth,  according  to  Singer. 

The  calculations  of  'this  section  make  no  pretense  to  accuracy,  but 
were  intended  primarily  as  a  vehicle  for  discussion  of  penetration  and 
shock  wave  models.  The  reader  is  also  referred  to>  recent  articles  by 
Axford  (1962)  and  Kellogg  (1962). 


E,  STABILITY 


Experimental  data  from  all  rocket  probes  so  far  have  indicated  a 
considerable  amount  of  turbulence  in  the  region  of  the  interface]  Ex¬ 
plorer  X  noted  close  inverse  correlation  vith  fluctuations  in  plasma 
density.  Even  if  the  Interface  is  stable,  the  fluctuations  in  the  SCR 
will  move  it  back  and  forth;  however,  it  has  long  been  believed,  on  the 
basis  of  other  geophysical  phenomena,  that  the  impact  of  sudden  fluctua¬ 
tions  generates  unstable  surface  waves  in  the  Interface.  Parker.  (1958) 
considered  the  magnetic  field  to  be  embedded  in  an  infinitely  conducting, 
incompressible,  inviscid  plsisma,  the  plane  surface  of  which  is  subjected 
to  a  wind  at  various  angles.  The  simple  application  of  linearized 
hydranagnetlc  theory  yielded  a  complex  dispersion  relation  for  waves 
in  the  surface,  and  convinced  Parker  that  instability  was  possible. 

Hxirley  (I96IC,  p.  42)  has  generalized  Parker's  approach  for  dif¬ 
ferent  orientations  of  magnetic  field  relative  to  the  impinging  wind 
and  finds  that  most  cases  of  interest  are,  in  fact,  stable.  However, 
in  the  noon  meridian  plane,  for  example,  instability  sets  in  as  one 
moves  toward  the  poles  past  a  certain  critical  latitude  9^  given  by 
(assuming  the  interface  approximately  circulsur  in  this  plane) 

tan  9^  =  2/6  ;  «  p/nm  ,  (104) 

where  p  is  the  density  of  the  plasma  below  the  interface  euid  nm  is  the 
density  of  the  SCR.  If  we  estimate  e  =  1  ,  then  9^  =  63.5° 
geomagnetic  latitude,  close  to  the  aurorsd.  zones.  The  curvature  of 
the  actual  field  would  add  to  the  stability. 

The  calculations  of  Parker  and  Hurley,  however  ideeillzed,  provide 
our  only  quantitative  basis  so  far  for  speculations  concerning  stability 
of  the  GMF-SCR  interface. 
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